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Abstract: 



We show that the Clifford group — the normaliser of the Weyl-Heisenberg 
group — can be represented by monomial phase-permutation matrices if and 
only if the dimension is a square number. This simplifies expressions for SIC 
vectors, and has other applications to SICs and to Mutually Unbiased Bases. 
Exact solutions for SICs in dimension 16 are presented for the first time. 
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1. Introduction 

The Weyl-Heisenberg group [1] first appeared in nineteenth century alge- 
braic geometry, and is at the roots of many things including harmonic analy- 
sis, theta functions, and — of course — quantum mechanics. Its automorphism 
group within the unitary group — the largest subgroup of the unitary group 
having the Weyl-Heisenberg group as a normal subgroup — appears in quan- 
tum information theory under the name of the Clifford group [2, 3]. 

The particular problem that motivated the present study is known as the 
SIC problem [4, 5]: in a complex Hilbert space of finite dimension N, find 
N 2 unit vectors \ipi) such that 

I^J =► |(^)| 2 = ^r- (1) 

Such a collection of equiangular vectors [6] is known as a SIC, which is short 
for a Symmetric Informationally Complete POVM, where POVM is short 
for Positive Operator Valued Measure. In physics a SIC represents a kind 
of fiducial measurement, of interest for quantum state tomography [5] and 
indeed for the very foundations of quantum mechanics [7, 8, 9]. Such mea- 
surements can be — and in low dimensions, have been [10, 11, 12] — realised 
in the laboratory, but — perhaps surprisingly — the theoretical SIC problem 
as stated is not easy to solve. At the moment exact solutions are known in 
dimensions 2-16, 19, 24, 35, and 48 (with 16 added here), while convincing 
numerical solutions are available in dimensions 2-67. These results are due to 
several authors; we refer to Scott and Grassl [13] for complete references, and 
just remark that beyond three dimensions the known SICs look appallingly 
complicated at first sight. 

Zauner's conjecture [4] states that in every dimension there is an orbit 
of the Weyl-Heisenberg group which forms a SIC, and moreover that every 
vector in such a SIC is left invariant by an element of the Clifford group 
of order three. Hence the problem of finding a SIC reduces to that of find- 
ing a suitable fiducial vector for the group to act on, and the second — very 
mysterious — part of the conjecture provides some guidance when one looks 
for such a fiducial vector. All available evidence supports Zauner's conjecture 
[13, 14]. 

The representation theory of the Weyl-Heisenberg group tells us that 
once one of its generators is given in diagonal form all its group elements 
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are represented by monomial phase-permutation matrices, that is by unitary 
matrices having only one non-zero entry per column, and per row [1]. Being 
products of a permutation matrix and a diagonal unitary, such matrices are 
also said to be of shift-and-multiply type. This property of the group can 
be traced back to the way that the representation is induced from that of 
its center. It is an important property shared by many, but not all, unitary 
operator bases of group type (or "nice error bases", as they are called in 
quantum information theory) [15]. 

In general the Clifford group is not represented by phase-permutation 
matrices. However, after a few preliminaries in section 2, we devote section 
3 to a representation of the Weyl-Heisenberg group which is special to the 
case when the dimension is a square, N = n 2 . It can be thought of as a 
finite dimensional analogue of the Zak basis [16], and is used in the theory 
of theta functions [17]. Our observation is that in this representation the 
entire Clifford group is given by phase-permutation matrices. In section 4 
we demonstrate that this remarkable feature is present if and only if the 
dimension is a square. In the remaining sections we explore some ways in 
which the phase-permutation basis can be useful. Section 5 illustrates how it 
is, in a way, aligned to SICs, while section 6 is devoted to exact solutions for 
SICs in 2 2 and 3 2 dimensions. In the former case they are trivial to obtain. 
In the latter they are not, but they look significantly better compared to how 
they look in the standard basis [13, 18]. In Section 7, we present new exact 
solutions for dimesion 4 2 ; previously these were known in numerical form 
only. Section 8 contains a remark on Mutually Unbiased Bases, and section 
9 summarises our conclusions. There are three appendices containing group 
theoretical theorems and the fiducial vectors for dimension 4 2 . 



2. Preliminaries 

We introduce the Weyl-Heisenberg group by writing down a defining repre- 
sentation. Choose a dimension N and assume that {|0), . . . , \N — 1)} is an 
ortho normal basis of C^. Define two phase factors w,r by 



3 



and two operators X, Z by 



X|u) = |u+1), Z\u)=u u \u), u E {0, . . . , JV — 1} , (3) 

where here and elsewhere the labels of the vectors are computed modulo JV. 
The matrix group generated by {r, X, Z} is the defining representation of the 
Weyl-Heisenberg group H(N). (It is known [1] that all irreducible represen- 
tations of H(N) in dimensions larger than one are unitarily equivalent to the 
defining one). Note that X and Z are represented by phase-permutation ma- 
trices, and the same is true for all the elements of the group since the product 
of two phase-permutation matrices is again a phase-permutation matrix. 

The Weyl-Heisenberg groups H(N) behave slightly differently depending 
on whether JV is even or odd. The underlying reason is that while one always 
has the relations 

u N = 1 , X N = Z N = 1 , (4) 
the order of r depends on the parity of JV: 

N _ J -1 JV even . . 

\ 1 JV odd. [b) 

As a consequence we will sometimes end up using arithmetic modulo N in 
the odd case and modulo 2JV in the even case. To unify the notation, we 
introduce the symbol 

/V _ / N N is odd 

\ 2N JV is even ' W 

Since ZX = uXZ, the introduction of the phase factor r may seem odd. 
If JV is odd, t is in fact a power of oj. The reason it is included in the 
even dimensional case can be traced back to the fact that there are group 
elements, such as XZ, that generate cyclic subgrops of order 2JV. For us 
it will be crucial that the Clifford group defined below acts on H(N) as we 
have defined it here. 

To analyze the structure of the Weyl-Heisenberg group, we define the 
group elements 

Dij = r ij X l Z j . (7) 
for i,j — 0, . . . , N — 1. One can then verify the central composition law 

DijDi m = r lj ~ im D i+ ij +m . (8) 
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Because the (non-scalar) generators X, Z of H(N) are just Diq and Dqi 
respectively, the composition law (8) implies that any element of H(N) is of 
the form r k Dij for suitable integers i, j, k. 

Note that the phase factor r 1 - 7 in (7) depends on ij modulo N, whereas 
X l Zi only depends on % and j modulo iV (by virtue of (4)). Hence the group 
law (8) says that H(N) modulo phase factors is isomorphic to Z^, where Z^v 
is the group of integers {0, . . . , N — 1} with addition modulo N. 

The situation encountered in the last paragraph will become a general 
theme below. Depending on our objective, we will take one of two points of 
view. Sometimes we will be concerned with concrete matrix representations 
of the groups involved. Then, we will be specific about all phase factors 
involved, and work with arithmetic modulo N. In other situations, however, 
a more abstract approach turns out to be beneficial. In these cases, we will 
factor out phases and work solely in terms of the discrete group I? N and its 
symmetry groups. 

To make the abstract approach more precise, let Z(N) be the center of 
the Heisenberg-Weyl group H(N). From the group law (8), it is evident that 
Z(N) = {r h l}fc =0 ,...,iv-i! elements are precisely the phase factors times 
the identity matrix. Our previous observation can now be phrased more 
succinctly as 

H(N)/Z(N) ~ Z 2 N . (9) 

Unless N is a prime, the integers modulo N form a ring but not a field. 
Therefore, strictly speaking, the "vectors" (i,j) G I? N are elements of a 
module rather than of a vector space. But we permit ourselves a slight abuse 
of terminology and speak of vectors in 1? N . 

We will be concerned with another group: the Clifford group. It consists 
of all unitary operators Uq normalising the Weyl-Heisenberg group, in the 
sense that for all i,j there are k' such that 

U G DijUl = r k 'D V j, . (10) 

Not all transformations i-> (i',f) are possible. The fact that the group 
law (8) involves addition of vectors in Z^ suggests that any such transfor- 
mation must be linear. Further, the fact that the symplectic inner product 
(Ij — im) modulo N of the vectors and (l,m) appears in (8) suggests 
that this inner product might be an invariant 

Ij — im = l'f — i'm mod N (11) 
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of the action of the Clifford group. 

These intuitions turn out to be true and yield an almost exhaustive un- 
derstanding of the Clifford group [14]. More precisely, recall that SL(2,N) 
is the group of linear transformations on Zjy leaving the symplectic inner 
product invariant. A Z^-valued matrix 

G=(Z ?1 (12) 



7 5 

is an element of SL{2, N) if and only if 

det G = a5 - £7 = 1 mod N . (13) 

In one direction, we have that for every G G SL(2,N), there is an element 
Uq in the Clifford group such that 

U G D l3 Ul = D Gi;y (14) 

A converse statement will be given below. The unitaries appearing in (14) 
are known explicitly (c.f. Ref. [14] for more details): if f3 is relatively prime 
to N — so that it has a multiplicative inverse — one finds 



U G = -±=e ld r^ Su2 - 2uv+av ^\u)(v\ , (15) 



N-l 

N 

u,v=0 

where e ld is an arbitrary phase. If (3 is not relatively prime to N we can use 
the decomposition 

a (3 \ ( -1 \ / -f + xa 5 + x/3 \ , . 

7 6 ) = { 1 x ) { -a -0 ) ' (16) 

where the integer x can always be chosen so that 5 + xf3 is relatively prime to 
N. We remark that in this representation symplectic matrices are represented 
by unitary phase-permutation matrices if and only if they are of the form 

G ={ a ia °>)- < 17 > 

They form a rather small subgroup. 
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Let us return to the more abstract point of view alluded to before. Since 
phase factors are left invariant by a unitary conjugation 

U G (r k l)Ul = r k t, (18) 

the Clifford group acts on H(N)/Z(N) ~ Z 2 N . It is easy to see that the action 
of the Clifford group on H(N)/Z(N) is precisely isomorphic to SL(2,N). 
This equivalence holds irrespective of whether N is even or odd and delivers 
the converse statement promised above. A proof is given in Appendix A. 

Symplectic matrices of order 3 are of special interest. If A^ > 3 it can be 
shown that a symplectic matrix is of order 3 if and only if its trace equals — 1 
mod N [14]. According to a precise form of Zauner's conjecture a SIC fiducial 
can always be chosen to be an eigenvector of the unitary Uz representing the 
symplectic matrix 




In the standard representation Uz is not given in monomial form (although 
monomial representations of other symplectic matrices of order 3 can be 
found for special values of A" [14]). 

The phase of Uz is chosen so that U\ is the identity operator. It is then 
possible, using Gauss sums, to deduce the dimension of the three eigenspaces 
of a Zauner unitary in arbitrary dimensions [4]. There are three eigenspaces 
£o>£i)£2 corresponding respectively to the eigenvalues 1, e 27 ™/ 3 , e 47 ™/ 3 . With 
an appropriate choice of phase one finds for their dimensions that 
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The numerical evidence [13] strongly suggests that the eigenspace S al- 
ways contains SIC fiducials, while il N = 3k or N = 3k + 1 the other 
two eigenspaces never do. The case N = 3k + 2 does not occur if N = n 2 . 

Finally we will be interested in the extended Clifford group, which in- 
cludes anti-unitary operators as well [14]. The determinants of the 2 by 2 
matrices are then allowed to take the values ±1. The extended Clifford group 
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divides the set of all SICs into orbits in a natural way: acting on a given SIC 
with an element of this group produces another SIC since we assume that 
the SIC itself is an orbit under the Weyl-Heisenberg group. 

3. Representation with phase-permutation matrices 

From now on we are in a Hilbert space of dimension N = n 2 , 

U N = U n ®U n . (20) 
We add one more phase factor to the ones we have defined already: 

uj = e N , r = —e N , a — e « . (21) 
Our key observation is that 

ZX = uXZ =4> [X n , Z n ] = . (22) 
Hence the Weyl-Heisenberg group admits an Abelian subgroup, 

S = (X n , Z n , rl) , (23) 

of maximal order NN. By diagonalising the elements of S, we define a new 
basis with basis vectors labelled |r, s), where r and s are integers modulo n. 
Indeed 

X n \r,s) = (T r \r,s) , Z n \r,s) =a s \r,s) . (24) 

Some phase choices are still to be made. We settled for a choice which implies 
that the generators of the group are represented by 

{|r, s + 1) if s + 1 7^ mod n 
(25) 
a r \r, 0) if s + 1 = mod n 

Z\r,s) = cu s \r - l,s) . (26) 

We refer to this representation as the phase-permutation representation. It 
treats the two generators X and Z in an even-handed way, so in a sense we 
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have gone "half-way" to the Fourier basis. The relation to the standard Weyl 
basis is given by 

1 n— 1 

\r, s) = —= V uT n * r |nt + s ) . (27) 

The matrix effecting this transformation is of the form F n <g) 1, where F n is 
the n by n Fourier matrix. 

The entire Weyl-Heisenberg group is represented by phase-permutation 
matrices, but this time more is true. In this representation the entire Clifford 
group is represented by phase-permutation matrices. The following armchair 
argument explains why: In all dimensions, the Clifford group permutes the 
various maximal Abelian subgroups of the Weyl-Heisenberg group. It also 
preserves the order of any group element. But if N = n 2 there is a unique 
maximal Abelian subgroup whose elements (modulo phases) have orders that 
are the divisors of n; namely the one that defines our basis. Therefore the 
Clifford group reorders the elements of this Abelian subgroup, and it follows 
that it reorders the basis vectors while possibly multiplying them with phases. 

The permutations involved are easy to deduce. Let G be a general sym- 
plectic matrix as given in eq. (12). Using its inverse we observe that 

U G X n U G = U G D n0 U G = Dsn-yn = T -' 5N X 5n Z^ n (28) 

U G Z n U G = U G D 0n U G = D_^ an = T -^ x -pn zan (2Q) 

As usual the case of odd n is simpler, since t n = 1 in this case. Let us 
therefore assume that n is odd to begin with. We see that 

X n U G \r, s) = U G X 5n Z^ n \r, s) = a 5r -~< s U G \r, s) (30) 

Z n U G \r, s) = U G X-P n Z an \r, s) = a-^ r+as U G \r, s) . (31) 

It follows that U G \r, s) is a common eigenvector of the diagonal operators X n 
and Z n , and indeed that 

n odd =4> U G \r, s) = e i6rs \5r - 7s, -0r + as) , (32) 
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where 6 rs is a phase to be determined. It again follows that an arbitrary 
symplectic unitary is represented by a phase-permutation matrix. 

The argument can be extended to the even dimensional case, but since 
we also need to calculate the phases 8 rs we will proceed a little differently in 
the general case. First we define 

{0 n is odd 
(33) 
| n is even . 

We may assume that the matrix element /3 is relatively prime to N, because 
if it is not we can fall back on the decomposition (16). Using the standard 
representation (15), and relation (27) between the two bases, it is straight- 
forward to show that 

K R-ltX a '2_o c , t ,lj_ nfC ,2\ 



(r',s'\U G \r,s) = I_ r ^W 2 -2-'W) x 

(34) 



n 2 



n-l 

v , ,nt(-r+p- 1 (-s'+as+ma)), ,nt' (r' +fi- 1 (8s' -s+m5)) 

X 7 UJ UJ j 

t,t'=0 

where m was defined in eq. (33). Performing the sums, and using the fact 
that a5 — f3 r y = 1 modulo n, we see that 



<r', s'\U G \r, s) = e^^-H^-^'W) ^ 



r = or — 7s — m - ^ ' 
s' = —fir + as + ma 



(35) 

and zero otherwise. If n is odd then m = 0, and we have reproduced eq. (32) 
but with the phases now included. If n is even we use modulo 2 arithmetic to 
polish the m- dependent term; note that (3 = 1 modulo 2 since /3 is relatively 
prime to n. Thus we arrive at our key result: 

Theorem 1: When the dimension is a square number N = n 2 , the Clifford 
group admits a representation by phase-permutation matrices. For the Weyl- 
Heisenberg subgroup it is given by eqs. (25). For an SL(2, N) element of the 
form (12), with (3 and N relatively prime and m and s' as above, it is 
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U G \r, s) = e< V 1 ( Ss ' 2 - 2ss ' +™ 2 ) \5r - 7 s + m 7 5, -fir + as + ma/3) . (36) 

Tie overall phase 8 remains undetermined. The case when f3 is not relatively 
prime to N can be recovered from eq. (16). 

Note that if G — G modulo n (using modulo n arithmetic for the matrix 
elements) then Uq and Uqi produce the same permutations of the basis ele- 
ments, that is to say they differ only by a diagonal unitary. 

The group element of most interest to us is Zauner's unitary, correspond- 
ing to the matrix (17). It is given explicitly by 

iir(N— 1) 2 I o 

U z \r,s) = e 12 r r +2rs \ - r - s - m,r) . (37) 

Here the overall phase 9 was chosen to ensure that U\ is the identity operator. 

A general element of the extended Clifford group is obtained by replac- 
ing SL(2,N) with the group ESL(2,N). The additional matrices E G 
ESL(2, N) can be written as a product 

E = GJ, GeSL(2,N), J=(^l _°i ) • ( 38 ) 

To the matrix J there corresponds an anti-unitary operator Uj whose action 
on the phase-permutation basis is given by 

Uj\r,s) = \-r,s) . (39) 

Hence, the extended Clifford group also acts through phase-permutation ma- 
trices on the basis vectors. 

4. Uniqueness proofs 

We have seen by construction that the Clifford group admits a representation 
by phase-permutation matrices if the dimension iV = n 2 is a square. We will 
now prove the converse, that this is possible only in square dimensions. Since 
the construction hinged on a special maximal Abelian subgroup of the Weyl- 
Heisenberg group — it transformed into itself under the Clifford group — we 
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begin with a theorem that shows that this is a necessary feature of any phase- 
permutation representation. It actually applies to a slightly more general 
situation, in which we consider a group Q which may be the entire Clifford 
group, but which may also be any subgroup of the Clifford group which 
includes the Heisenberg group H(N) subgroup. 

An Abelian subgroup S of H(N) is called maximal if no element g G 
H(N) \ S commutes with everything in S. Equivalently, S is maximal if it 
has the maximal possible order NN. 

Theorem 2: The following two statements are equivalent: 

1. There exists a phase-permutation representation of Q on C N which is 
irreducible when restricted to H(N). 

2. There exists a maximal Abelian subgroup of H(N) which is stabilised 
byG. 

Proof. First we establish 2 =>• 1. Let S be a maximal Abelian subgroup of 
H(N). Choose a basis in which all elements of S are simultaneously diagonal. 
If we select one representative from every coset in S/Z(N), their diagonals 
define A" orthogonal vectors (because H(N)/Z(N) defines a unitary operator 
basis). From this one concludes that a maximal Abelian subgroup defines 
a joint eigenbasis which is unique up to permutations and rephasings. It is 
then obvious that 2 1. This was used in section 3. 

To prove that 1 =>- 2, denote the phase-permutation basis by {\e a )}^ =1 . 
By assumption H(N) simply permutes the corresponding rays (vectors up 
to phase), and acts transitively on them. Let S be the subgroup of the 
Heisenberg group leaving the particular projector |ei)(ei| invariant. Since 
the orbit of H(N) acting on the set of projectors has size A" we know that 
\H(N)\/\S\ = N, that is to say that 

\H(N)\ N 2 N 

where N is the cardinality of the center. Modulo phases, S must thus be 
a subgroup of H(N) of order N. By definition, all elements of S have the 
common eigenvector \ei). But it is a direct consequence of the commutation 
relations that two elements of the Heisenberg group have a common eigen- 
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vector if and only if they commute. Therefore S is an abelian subgroup of 
H(N) and has order NN. Hence S is maximally Abelian. 

Now H(N) acts monomially on the joint eigenvectors of any maximal 
abelian subgroup. It also acts transitively, by irreducibility. Thus the orbit 
of |ei) under H(N) consists precisely of the joint eigenvectors of S (up to 
phases). But we showed before that the orbit coincides with {\e a )}^ =l . □ 

We now focus on the case where Q is the full Clifford group. Are there 
cases beyond square dimensions where there is a maximally Abelian subgroup 
S of H(N) stabilized by the Clifford group? By Section 2, the orbit of an 
element of H(N)/Z(N) under the action of the Clifford group corresponds 
to the orbit of a vector, v G 7? N , under the action of SL(2, N). The following 
characterization of these orbits is implied by Lemma 27 of [19]. We re-prove 
it here to make the presentation self-contained. 

Recall that the order of a vector v G I? N is the least integer k > 1 such 
that kv = 0, where the triple bar equality sign denotes equality modulo N. 

Lemma 1: The orbits of the action of SL(2, N) on 7? N are the sets 



of vectors of constant order. 

Proof: We first prove that S G SL(2, N) cannot change the order of a vector 
v. Indeed, if kv = 0, then 



so that ordSt> < ordi>. Replacing S by S -1 , we see that ordSt> > ordv. 
Hence equality must hold as claimed. 

The comparatively difficult part is to show that SL(2, N) acts transitively 
on the sets O k . Let v G Z 2 N be of order k. By definition we thus have that 



O k = {v G Z 2 N | ordv = k} 



(41) 



k(Sv) = S(kv) = SO = 



(42) 



N I ( kv i) ^~r\ v h 



(43) 



where Vi,i = 1,2 are the components of v. Therefore, the vector 



1 



(44) 



N/k 



v 
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is well-defined as an element of T? N . One checks that ordf' = N. 

We go on to show that there is a symplectic matrix S G SL(2, N) whose 
first column equals v' . That is the case if there are integers x,y such that 

1 = det ( v ) 

By Bezout's identity, there are integers a, b such that 

v[a + v' 2 b = g, (46) 

where g = gcd(v[, v' 2 ). It must be the case that g and N are co-prime, for 
otherwise (N/g) would be an integer smaller than N such that (N/g)v' = 
0, which would contradict the fact that ordv' = N. Thus there exists a 
multiplicative inverse g^ 1 of g modulo N. Hence 

V = g^a, x = -g^b (47) 

provides a solution to (45). 

Finally, let w be another vector of order k. Let S v be a symplectic matrix 
with first column equal to v', let S w be a symplectic matrix with first column 
equal to w'. Then 

SwS-W = w'=> S w SZ x v = w. (48) 

Thus any two elements of Ok can be mapped onto each other by means of 
an element of SL(2,N). □ 

The preceding lemma allows us to decide in which dimensions there is a 
monomial representation of the Clifford group just by counting orbit sizes. 
It seems simpler to do that in prime-power dimensions. 

Lemma 2: Let N — pf . . . p q k k he the decomposition of the dimension into 
powers of distinct primes. 

There is an order-N subgroup of 1? N which is stabilized by SL(2, N) if 
and only if the same is true for all dimensions iVj = pf, for i = 1, . . . , k. 

The statement follows from the more general fact that the Weyl-Heisenberg 
group and the Clifford group factor into direct products for composite N. It 
is proven in Appendix B. We are ready to conclude this section: 



VnX. 



(45) 
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Theorem 3: There exists a monomial representation of the Clifford group 
which contains the Weyl-Heisenberg group as an irreducible subgroup if and 
only if the dimension N = n 2 is a square. 

Proof: Using the notions of Lemma 2, let iVj = pf. Let V C 1? N . be a 
non-trivial subgroup which is invariant under the action of SL(2, N). 

Let k = max{ordf \v G V} be the largest order of any element in V. 
By Lagrange's Theorem, k is of the form k = p\ for < I < qi. Because 
v G V =>- piV G V, there is also an element of order p l ~ l in V, and, indeed, 
any power of pi up to the /th appears as the order of some element in V. 

By Lemma 1 and the assumption that V be SL-invariant, we find that 

V = {veI? Ni \ordv<p\}=pf- l Zl, (49) 

Hence |V| = pf, which is equal to TVj if and only if / = qi/2. That is possible 
if and only if q^ is even, which implies the claim. □ 

In the remaining sections we turn our attention to applications of the 
phase-permutation basis. We return to our original motivation of SICs. First 
we discuss a general property of SICs on the phase-permutation basis, and 
then in sections 6 and 7 we use this basis to construct SICs in dimensions 
2 2 , 3 2 , and 4 2 . The sixteen dimensional case has so far not been solved in 
the Weyl basis. Finally, in section 8 we consider sets of Mutually Unbiased 
Bases (MUB) in the phase-permutation basis. 

5. Images of SICs in the probability simplex 

For the first application, recall that a quantum state together with an or- 
thonormal basis gives rise to a probability vector, with the basis vectors 
themselves forming a simplex. In physics this is connected to quantum mea- 
surements, in mathematics it is an example of the moment map. Anyway, if 
the quantum state is a pure state represented by the Hilbert space vector 

(z , ^_x) T = (V^, Vple tfl \ . . . , VmTie^f , (50) 

then its image with respect to the same orthonormal basis is a probability 
vector with components pi. The probability vector gives the barycentric 
coordinates of a point within the probability simplex associated to the basis 
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chosen. We assumed that the vector is a unit vector, which indeed implies 
that 



If we have a set of such pure states connected by a group which is represented 
by phase permutation matrices, then their probability vectors will be related 
by permutations of the coordinate axes. It follows that they lie on a sphere 
centered around the midpoint of the probability simplex. 

More can be said if we deal with the N 2 unit vectors forming an orbit 
under the Weyl-Heisenberg group and obeying the SIC condition. Let the 
components of a SIC fiducial vector be given by eq. (50). In the standard 
basis, where the subgroup generated by Z is diagonalised, it is known that 



The length of the probability vector, that is the radius of the sphere in which 
it is inscribed, is thereby determined. The remaining equations also admit 
an interesting geometrical interpretation [22]. Since the matrix representing 
Z is diagonal, the N states Z r \ip F ) all project to the same probability vector. 
Moreover the ith component of the vector X x \%Pf) is Pi- X , so from these 
equations we can read off not only the length of the probability vectors but 
also their mutual angles. The conclusion turns out to be that the image 
of the N 2 SIC vectors is itself a regular simplex with N vertices. Once this 
geometrical interpretation is available it is unsurprising that the N equations 
(52) are redundant and do not by themselves determine all the coefficients 
Pi (unless N = 2). But they are still helpful. 

The reason why the iV 2 vectors in the SIC give rise to only iV images 
in the projection is that the images form an orbit under the subgroup that 
is complementary to the diagonalised subgroup. It is clear that something 
similar should happen in the phase-permutation basis, where again there is 
an Abelian subgroup of order N that does not move the projected points. 
To show this we denote the components of the SIC fiducial \i^f) by z TS = 
Jp^lz 1 ^ 3 , and find 





(51) 



[20, 21] 




(52) 
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n-1 

(i, \X nU Z nV \^)= Y,Prs° rU+SV 
r,s=0 
n—l n—l 

\(4> \X nU Z nV \4> )\ 2 =J2Y1 PrsPr's>V {r - r ' )U+iS - S ' )V ■ (53) 
r,s=0 r',s'=0 

The absolute values on the left hand side are known from the condition 
defining a SIC. Using this, and summing over the integers u and v, we find 

S^ = iv( 1+ ivTTj = ivTT- (54) 

More generally we can take a Fourier transform of eq. (53). This gives N — l 
additional equations for the absolute values, 

n—l j 

^ PrsPr+x,s+y = ~ ^ • (55) 
r,s=0 

Here x, y are integers modulo n, not both zero. This is analogous to what 
happens in the standard basis, and the geometrical interpretation is the same: 
when the SIC is projected to the basis simplex we see a regular simplex 
centered at the origin with just N vertices. Its orientation differs from the 
projection to the standard basis (see Fig 1). 

6. SIC fiducials for N = 2 2 and 3 2 

If we use the phase-permutation representation it becomes very easy to find 
the SICs in N = 2 2 dimensions. In fact the absolute values entering the 
fiducial are determined by Zauner's conjecture, normalisation, and eq. (54). 
Alternatively one can simply solve the equations defining a Weyl-Heisenberg 
covariant SIC. Before doing so it is convenient to rephase the basis through 

(|0,0), |0,1), |1,0), |1,1» (r- 2 |0,0>, t- 7 |0,1>, t- 5 |1,0>, |1,1». (56) 
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Figure 1: The images of an N = 4 SIC, in the two bases we discuss. The 
phase-permutation basis (the large tetrahedron on the right) is aligned to the SIC 
(the small tetrahedron) in a way that the standard basis (the large tetrahedon 
on the left) is not. The corresponding pictures for N = 9 (not shown) are eight 
dimensional, and would show that the faces of the eigenvalue simplex for the 
phase-permutation basis are nicely aligned with the image of the SIC. 



This gives the representation 
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(57) 



There are altogether 4 4 = 256 possible solutions for a SIC fiducial, giving 
rise to 16 SICs altogether. The solutions are 



\ 



/ x 

i si 

\ i Ul ) 



( i* 2 \ 

x 

\ i U2 J 



( i* 3 \ 



x 
m 



( i s * \ 



\ i U3 ) 



,t 4 



(58) 



V x J 



where Si,ti,Ui are integers from to 3, 



x 



2 + V5 



(59) 



and an overall normalisation has been ignored. This simple form of the SIC 
vectors has been found before by Zauner [23], who arrived at it by casting 
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the Zauner unitary into phase-permutation form, and from an alternative 
point of view by Belovs [24] . Each of the vectors is left invariant by a Zauner 
unitary of order 3, and the 16 SICs form a single orbit under the Clifford 
group. The result of the numerical searches [5, 13] is thus fully confirmed. 
The images in the probability simplex (see section 5, and Fig. 1) coincide for 
all the SICs, and they are nicely oriented. The structure of the N = 4 SICs, 
and their entanglement properties, were studied in detail recently [25]. 

The case of N = 3 2 is much harder. Lest our readers be disappointed by 
this, we strongly recommend they begin by looking at the answer found — 
by means of a Magma calculation — in the standard representation [13, 18]. 
Afterwards our result will come as a pleasant surprise. 

With our canonical choice of Uz, Zauner's conjecture implies that the 
SIC fiducial takes the special form 

|V>) = -,W|1,1) -z 2 u\2,2) + 

(60) 

+^ 3 (w 6 |0,2) + |1,0) +c/|2,l)) + z 4 (u 6 \0,l) + |2,0) +w 5 |l,2)) . 

We included some convenient phase factors. We exploit the arbitrariness in 
the overall phase to write the Zj in the form 

Zl = VP~ie m) Z2 = Vp~2e~ m> z 3 = v fae i » 3 z 4 = ^ 4 e^ (61) 

where it is assumed that — n/2 < fi < ir/2. 

The necessary and sufficient condition for a normalized vector \ip) to be 
a fiducial vector is that 

im-^>i 2 = ^ (62) 

for all j, k not both zero — a total of 80 equations. However the Zauner sym- 
metry means that if eq. (62) is satisfied for the vector (j, k) it is automatically 
satisfied for the (up to) two other vectors obtained by acting with the Za- 
uner matrix. Also, if it is satisfied for (j, k) it is automatically satisfied for 
(— j, —k). Consequently we can reduce the 80 equations to 15, which it will 
be convenient to group as follows 
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1. 2 group 1 equations 

l<^Po^>| 2 = |^|I>s,6^>| 2 = ^ (63) 
Assuming normalization these equations are equivalent to eqs. (55). 

2. 3 group 2 equations 

MD 3j+1:6j+2 m 2 = ± (64) 

with j = 0,1,2. 

3. 9 group 3 equations 

MDw^Mf 1 = ± (65) 

with j,k = 0,1,2. 

Together with normalization this gives us 16 conditions on the 7 real parame- 
ters in eq. (61), so there is still a high degree of redundancy in the equations. 
The nature of the dependencies will become clear in the course of solving 
them. 

We begin by considering the normalization condition and group 1 equa- 
tions. As can be seen from eqs. (54) and (55) (which are equivalent to them) 
they only involve the absolute values. Explicitly: 

P1+P2 + 3p 3 + 3p 4 =1 

P 2 i + p\ - P1P2 = jo • (66) 

3pg + 3 P4 + 3p 3 P4 -Ps-Pi = ~To 
These equations are not hard to solve. Setting 

pi = ai + h , p 2 = ai - 61 , p 3 = a 3 + 63 , p A = a 3 - b 3 (67) 
diagonalizes them. It is then readily deduced 
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a 3 = i(l-2ax) , b\ = ^(l-10a?) , b 2 3 = -L (-1 + 20^ -60a?) . (68) 

To fix the free parameter in the expressions just derived we need to con- 
sider the group 2 equations. This will also give us the phase e lMo . It is 
convenient to write the equations in the form 



r,s=0 

with j = 0,1,2 and where 



E ^~% r e: = 1 (69) 



e — z l z 2 i e l 



= iV3co 2 p 3 , e 2 = — iv^w 7 ^ . (70) 



Inverting the Fourier transform we see that eqs. (69) are equivalent to the 2 
conditions 



\„ |2 i U |2 | i |2 1 

|e | + |ei| + |e 2 | = — 

(71) 

e e l e l e 2 ~l~ e 2 e = 

The first of these is a consequence of the group 1 equations. Solving the 
second for e = z x z\ gives 

* _ e 2 2 el - efe^ _ V3(a 2 - b 2 )(V3b 3 - ia 3 ) 

12 " h| 2 - |e 2 | 2 - 4a 3 6 3 ' ' ( j 

Taking the square of the absolute value on both sides and using eqs. (68) we 
find 



(1 - 40a! + 40a 2 )(-ll + 100a! - 120a 2 - 800a? + 1600a?) = . (73) 

Solving this equation and taking account of the requirement that bi, b 3 both 
be real we deduce 
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ai = ^5 - s 5\/3 + s 3VE + 



60 



-18 - s 7V3 + s 6VE + 5VTbj 



(74) 



where Sq, si and s 2 are arbitrary signs. This fixes the absolute values. Note 
that the only choice of sign that affects the set of absolute values is s , which 
suggests — correctly — that s labels two different orbits of the Clifford group. 

To determine the phase e w , we substitute these expressions into eq. (72) 
and simplify. We obtain 



e 2ifM> = i ^2 (6 + s V3-Vlt) - ^2 (2 - s V3 + VTE^j (75) 
Taking account of the assumption that — ir/2 < pL < ir/2 we deduce 



e™ = )Jl + c -i8i ) J\-c (76) 

where 

00 = ^2(6 + ^3-^) . (77) 

Note that the numbers given in these expressions as nested square roots can 
be constructed with ruler and compass, so the ancient Greeks might have 
approved — especially since the 9th root of unity cannot be so constructed. 

To calculate the remaining two phases we turn to the group 3 equations. 
It is convenient to write the equations in the form 

E v j{r ~ s) e kr el s = ± (78) 

r,s=0 
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z 3 z 4 



with j, k = 0, 1, 2 and where 

e fc0 =(l + 2(-l) fe cos^) 
e kl = - (r 6k ~ 5 zl + t'^-V z* 2 ) z 3 

e k2 = - (T 6k - 5 Zl + t-^-^Zz) z\ . 

Inverting the Fourier transform in eqs. (78) we see that the nine group 3 
equations are actually equivalent to the six equations 

|e fc0 | 2 + |e fe i| 2 + |e fc2 | 2 = (79) 

e fco4i + e ki e *k2 + e k2 e lo = ( 80 ) 

with k — 0, 1, 2. Using eq. (72) and some elementary trigonometry one finds 
that for all three values of k eq. (79) is equivalent to the single condition 

3(al-bl)+4 ai a 3 = ^ (81) 

which is an immediate consequence of the group 1 equations. We are thus 
left with the three eqs. (80). It will be convenient to write them in the form 

h = h = h = (82) 

where 

2 

fj = ° 3k ( e fco4i + e kl e* k2 + e k2 e* k0 ) . (83) 

fc=0 

Writing the expressions out in full we find that / 2 = r 2 f\. So the nine 
equations with which we started reduce to just the two equations fo — fi — 0. 
It is readily confirmed that these are equivalent to 
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3iM3 _ \ ( ^1 + Z l Z 2 + Z 2 |_ 1 ( Z l Z 1 Z 2 + Z 2, 



gOifj,3 I ± f _|_ 



(84) 

1 / z\ + + % {z\ - z x z 2 + zf) 



2^pl\ zl + z* y/3(zt-zZ) 



The quantities on the right hand sides are all known so these formulae give 
explicit expressions for the two remaining phases. Simplifying them and 
taking the cube roots we find 




(85) 



e »M4 _ a m 4 




where 



ci = ^V9- s W3 + 503^-2^15 
8 

SiSq 

24~ 



(86) 



c 2 = 



\J 15(— 19 + so^v^ - s 9y/5 + 6-/15) . 



Here m 3 and m 4 can take the values 0, 1,2. The entire solution is given in 
terms of radicals, as expected (but not understood!). 

Altogether there are 2 3 • 3 2 = 72 fiducial vectors, splitting into 2 different 
orbits of the extended Clifford group labelled by so — ±1. The solution with 
so — s\ = S2 = 1TL3 = = 1 is the fiducial 9a as labelled by Scott and 
Grassl [13], while switching the sign of (only) s leads to their fiducial 9b. 
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Figure 2: A change of basis for dimension N = 16 resulting in a Zauner 
matrix which is a permutation, where r = — exp(7r«/16). 

7. SIC fiducials for N = 4 2 

The first dimension for which the approach helped in finding a new solution 
is iV = 16 = 4 2 . We were unable to obtain a solution by hand as we did for 
N = 9. However, we were able to obtain a solution using Magma. 

We use a basis such that both X 4 and Z 4 are diagonal. The change of 
basis is given by the matrix T in Figure 2, and the Weyl-Heisenberg group 
in this basis are shown in Figure 3. 

In this basis, the Zauner matrix is a permutation matrix. Hence a fiducial 
vector is of the form 

|V) = x (|0) + |2> + |6» + Xl (\l) + |9> + |10» + x 3 (|3) + 1 14) + |15» 

(87) 

+ x 4 |4) + x 5 (|5) + |11) + |12)) + x 7 (\7) + |8) + |13)). 

In order to solve the equations for these six complex variables Xi, one of which 
can assumed to be real, we followed the approach described in [18]. Com- 
puting a Grobner basis modulo a single 23-bit prime using Magma [26] took 
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Figure 3: Generators X and Z of the Weyl-Heisenberg group in dimension 
N = 4 2 in the adapted basis. Here r = — exp(rri/16) denotes a primitive 
32nd root of unity. 
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about three days and required about 30 GB of memory. The polynomials in 
a Grobner basis with respect to so-called grevlex order have coefficients with 
90 digits in the numerators and denominators. Changing to lexicographic 
order which is used to solve the equations, the coefficients grow to some 900 
digits. Nonetheless, we succeed to obtain a less complex representation of a 
solution. The solutions are given in a number field 

(88) 



K = Q(v / 2, Vl3, Wf, r 2 , r 3 , t u t 2 , t 3 , t 4 , y/=T), 
of degree 1024, where 



r 2 
ti 

4-2 



= \ V221- 11, 



?~3 



\/l5 + VV7, 



^15 + (4 - vTr> 3 - sVu, (89) 

(((3 - 5^)^ + (39^ - 65))r 3 + ((16>/i7 - 72)^ + 936))*i 
- 208^ + 2288, 

t 3 = ^2-V2, U = V^+h- 

Note that we do not explicitly use a 32nd root of unity, it can be expressed 
as 

u 32 = ^((V2(l-t 3 )-l)U-UV^l). (90) 

The Galois group of K is isomorphic to C 8 x ((C 2 x C 2 x C 16 ) x C 2 ), and 
K is an Abelian extension of Q(v / 22l). The coefficients of a non- normalized 
fiducial vector of eq. (87) are as follows: 

x = - ^>/l3r 3 tit 2 (91) 
.ti = (^((21v^+22v^+ 16v^ + 5v^+5v^M + 4v / 22T+v / 442 + 74)r 2 r 3 

+ (-77V2 - 26\/l3 - 18VT7 - 33V26 - 19V34 + 2V221 - 7V442 + 42)r 2 
+ (-45V2 + 30VT3 - l0Vl7 + 15V26 + 5V34 + 10V221 + 5V442 - 30)r 3 
+ (30Vl3 + 30Vl7 + 10V22T - 70))ti* 4 

+ ((3V2 + 3y/l3 + 9VV7 - V26 + 7V34 + 11 V221 + 3V442 + 121)r 2 r 3 (92) 

+ (-82V2 - 88\/l3 - 24V17 - 74V26 - 2VM + 24V221 - 2\/442 + 264)r 2 

+ (175V2 + 80Vl3 - 20Vl7 + 75V26 - 25V34 - 5V442 + 380)r 3 

+ (200V2 + 220V13 - 300VT7 + 160V26 - 160V34 - 20V221 - 40V442 + 180))t 4 ) ' 
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£3 



+ 
+ 
+ 
+ 



+ 
+ 
+ 
+ 
+ 

+ 
+ 
+ 

+ 
+ 
+ 
+ 

Xi = ((■ 
+ 
+ 
+ 



(-10V2- 15-/L3- 15V17- 6V26- 8V34- \/221- 21)r 2 r 3 

55\/2 + 16-/L3 + 28-/L7 + 7V26 + 21 V34 + 8\/22T + 5^442 + 108)r 2 

70V2 + 10V26 + 80)r 3 

10^2 - 130v/l3 - 250^17 - 70\/26 - 70\/34 - 30V221 - 10V442 - 630)tii 4 
(10\/2 - 51\/l3 - 33V17 - 24^26 - 22\/34 + V221 - 29)r 2 r 3 
320^ - 4\/l3 + 28\/l7 + 8V26 + 4V34 + 44v / 221 + 20\/442 + 524)r 2 
265^2 - 30Vl3 - 50\/l7 - 15V26 - 35^34 + 10^221 + 5\/442 + 310)r 3 
260A/2 - 200^13 - 560\/l7 - 100^26 - 260^34 + 20^442 - 600))i 4 
-21\/2 - 22-/L3 - 16\/l7 - 5^26 - 5%/34 - 4\/221 - VU2 - 74)r 2 r 3 

77 + 26-/L3 + 18-y/l7 + 33V26 + 19V34 - 2\/221 + 7\/442 - 42)r 2 
-45V2 + 30V13 - 10\/l7 + 15V26 + 5\/34 + lO-y/221 + 5\/442 - 30)r 3 
30Vl3 + 30\/l7 + 10V221 - 70))ti* 4 

(-3\/2 - 3V13 - 9VT7 + V26 - 7V34 - 11V221 - 3\/442 - 121)r 2 r 3 

82V2 + 88\/T3 + 24\/l7 + 74V26 + 2\/34 - 24\/22T + 2\/442 - 264)r 2 

175\/2 + 80v/l3 - 20Vl7 + 75V26 - 25^34 - 5\/442 + 380)r 3 

200V2 + 220V13 - 300V17 + 160V^6 - 160^ - 20V221 - 40V442 + 180))t 4 ) 

(10\/2 + 15-/L3 + 15\/l7 + 6V26 + 8V34 + V22l + 21)r 2 r 3 (93) 
-55\/2 - 16\/l3 - 28\/l7 - 7^26 - 21^34 - 8^221 - bVU2 - 108)r 2 
70V2 + 10V26 + 80)r 3 

10^2 - 130V13 - 250\/l7 - 70^26 - 70V34 - 30^221 - lO-y/442 - 630)tii 4 
(-10V2 + 51V13 + 33Vl7 + 24V26 + 22V34 - V221 + 29)r 2 r 3 
-320V2 + 4VT3 - 28V17 - 8V26 - 4VM - 44^221 - 20^442 - 524)r 2 
265^ - 30Vl3 - 50VV7 - 15\/26 - 35^34 + lOv/221 + bVU2 + 310)r 3 
260V2 - 200\/l3 - 560\/l7 - 100\/26 - 260^34 + 20^442 - 600))i 4 

— Vl3- -Vl7- -x/221- -W 3 
26 2 26 2' 6 



lO-y/2- 



20 
13 



26 



10 ^ 
13" 



/ 442)r 2 )t 1 t 2 V-l 



(94) 



(|VT3 + i^+ 1^221 + i)r 2 , 3 
10V2+ §V26+ ^V442)r 2 )iit 2 
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+ 
+ 
+ 
+ 

(( 

+ 

+ 
+ 
+ 
+ 
+ 
+ 



75 V2 - 4^13 - 2VT7 - ? V26 - ^Vsl 



2 2-" 2— fVi^+lO)^ 

-22V2 - 24VT3 - 12VT7 + 14V26 + 2V34 - 4\/221 - 2V442 - 24)r 2 
15V2 - 5V13 - 5\/l7 + 25V26 - 5\/34 - 5V221 + 5\/442 + 35)r 3 
270^2 + 60^13 + 180VT7 + 10V26 + 70V34 + 20V221 + 10V442 + 620))ii 
(-85^2 - 28\/l3 - 4\/T7 - 3\/26 + V34 + 4^221 - 5\/442 - 36)r 2 r 3 
-190^2 - 86V13 + 22\/l7 + 34\/26 + 22V34 + 22\/221 - 10V442 + 122)r 2 
300V2 - 60\/l3 + 40VT7 + 40V26 - 20^34 - 220)r 3 
650\/2 - 60\/l3 + 460^17 + 110\/26 - 130\/34 + 60\/221 + 10V442 + 660 

(li\/2 + 23\/l3+ 19VT7- ^V26- ^V34 + 3 \/22T+ i\/442 + 63) r 2 r 3 (95) 

152^2 - 20^17 + 28\/26 + 24V34 + 4\/221 + 12V442 + 64)r 2 
-70V2 - 5\/l3 + 15VT7 + 10V34 - 5^221 + 55)r 3 

350V2 - 100^13 - 100\/T7 + 50V26 + 110V34 - 20V221 + 10\/442 + 60))ii 
43V2 + 28\/l3 + 24\/l7 - 23V26 - 19 VM + 8\/221 + 3V442 + 108)r 2 r 3 
476V2 - 22\/l3 - 26^17 + 28^26 + 4\/34 + 6\/22T + 36^442 + 26)r 2 
-170V2 - 20\/l3 - 40\/l7 + 50V26 + 10V34 - 10\/442 + 60)r 3 
410\/2 - 160VT3 - 120V17 + 150^26 + 270\/34 - 30VU2 + 280 

-li\/2 - 23\/l3 - 19VT7+ ^V26+ ^\/34- 3\/22T- i V442 - 63)r 2 r 3 

-152V2 + 20Vl7-28%/26- 24V34-4V22T- 12^442 - 64)r 2 
-70V2 - 5\/T3 + 15VT7 + 10V34 - 5^221 + 55)r 3 

350\/2 - 100^13 - 100\/l7 + 50\/26 + 110\/34 - 20^221 + 10^442 + 60))ii 
(-43^2 - 28\/l3 - 24\/T7 + 23V26 + 19\/34 - 8\/221 - 3\/442 - 108)r 2 r 3 
-476 V2 + 22V13 + 26V17 - 28^26 - 4\/M - 6>/221 - 36^442 - 26)r 2 (96) 
-170V2 - 20Vn - AOVU + 50^26 + 10V34 - 10n/442 + 60)r 3 

410V2 - 160\/l3 - 120\/l7 + 150V26 + 270V34 - 30V442 + 280))) V^T 

75 25 15 5 

(-y\/2-4VT3- 2\/l7- yV26- y\/34- -V442 + 10)r 2 r 3 



-22^ - 24\/l3 - 12\/l7 + 14\/26 + 2^34 - 4^221 - 2\/442 - 24)r 2 
- 15 V2 + 5\/l3 + 5V17 - 25\/26 + 5 V34 + 5V221 - 5 V442 - 35)r 3 
- 270V2 - 60Vl3 - 180Vl7 - 10V26 - 70V34 - 20^221 - 10V442 - 620)*i 
+ (-85V2 - 28VT3 - 4-/L7 - 3\/26 + V34 + 4\/22T - 5\/442 - 36)r 2 r 3 
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+ (-190V2 - 86\/l3 + 22\/V7 + 34V26 + 22^34 + 22^221 - 10V442 + 122)r 2 

+ (-300V2 + 60Vn - AOVU - 40V26 + 20V34 + 220)r 3 

- 650\/2 + 60Vl3 - 460VT7 - 110\/26 + 130\/34 - 60%/22T - 10V442 - 660 

Fiducial vectors for SICs with respect to the standard representation of the 
Weyl-Heisenberg group are given in Appendix C. It turns out that the two 
numerical solutions with orbits labelled 16a and 166 in [13] are related by 
the Galois automorphism of K induced by simultaneously changing the signs 
of and y/Tf. 

8. A remark on Mutually Unbiased Bases 

Although our main emphasis has been on SICs, we note that square dimen- 
sions are special also in the Mutually Unbiased Bases (MUB) problem. In 
dimensions N = n 2 , Wocjan and Beth have shown that one can construct 
sets of Mutually Unbiased Bases using n by n Latin squares [27]. In this 
section, we examine a variant of the MUB problem making use of the phase- 
permutation basis. 

We would like to find all vectors, and all bases, unbiased with respect 
to the two eigenbases defined by two complementary cyclic subgroups of the 
Weyl-Heisenberg group. This problem has been studied in connection with 
the MUB existence problem in dimension six [28] and in another guise is also 
known as the cyclic A- roots problem [29]. The two bases can be taken to be 
the standard basis \u)q, the eigenbasis of the subgroup generated by Z, and 
the Fourier basis |w)oo, which is the eigenbasis of the subgroup generated by 
X. (The labels and 00 do have a logical explanation [21].) We are looking 
for all vectors such that 



for all A values of u. The answer is known for N < 9, and it is also known 
that such vectors always belong to a complete basis unbiased with respect 
to the two bases we start out with. Here we just want to report what this 
problem looks like when N = n 2 and the phase-permutation basis is used. 




(97) 
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The two eigenbases are now given by 



n—l 

a + nb) = -=V a- br \r,a) , Z\a + nb) = u a+nb \a + nb) (98) 
\Jn ^— ' 



r=0 

n—l 



1 '* " 

la + nfo)^ = ^^cr 6r ur a >,r) , X|a + n6)oo = w a+nb |a + n6)oo • (99) 
We look for bases unbiased with respect to this pair, of the form 

j n—l 

|a + n6) fc = — ^2u k (r,a,b)\r, X k (r,a)) . (100) 



n 

r=0 



Here Wfc(r, a, 6) is a phase factor, and \ k is a map from Z„ x Z n to Z n , where 
Z„ is the set of integers modulo n. 

Unbiasedness with respect to the standard basis holds if and only if the 
function \ k (r,a) is injective for fixed a. To see this, we argue by contra- 
diction. Suppose Afc(r, a) is not injective for some fixed a, then there exists 
an integer x G Z n such that A^(r, a) 7^ x for all r. Now consider the inner 
product 



n—l 



(x + nb'\o a + nb) k = -^2q b ' r u k (r,a,b)S ZjXk (r,a) (101) 

r=0 

= 0, (102) 



since by assumption, S Xt x k ( r ,a) — for all r. Hence, if is not injective there 
is at least one vector from the standard basis that is orthogonal to the new 
basis. There are is no way to choose the phases u k to make it MU to every 
vector in the standard basis. 

Similarly, unbiasedness with respect to the Fourier basis holds if and only 
if the function is injective for fixed r. This means that the function X k (r, a) 
defines a Latin square, an n by n array such that each row and column 
contain the symbols from an n-letter alphabet exactly once. 

As an example, let n = p be a prime number. Then the Weyl-Heisenberg 
group contains p+1 cyclic subgroups altogether. Two of them were accounted 
for from the start, and the remaining p — 1 examples give rise to the choice 

\ k (r,a) = a + kr , k e {1, 2, . . . ,p - 1} . (103) 
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By inspection one finds that these p — 1 bases are maximally entangled. In 
addition, the Latin squares that define the bases are Mutually Orthogonal, 
which ensures that we have a collection of p + 1 Mutually Unbiased Bases. 

The fact that x — 1 Mutually Orthogonal Latin squares of order n give 
rise to x + 1 Mutually Unbiased Bases — regardless of whether they originate 
from group theory or not — was first observed by Wocjan and Beth [27]. The 
only new observation here is that Latin squares appear naturally in the cyclic 
iV-roots problem when the phase-permutation basis is used. 

To illustrate the idea, consider dimension N = 2 2 . The two eigenbases 
are given by 



and 



|0>o = 
|l)o = 

|0)oo = 
IDoo = 



^(10,0) + 



1 

7! 



(|o,i) + 



= ^(|o,o> + 
7! (|1 ' 0) - ? 



The Latin square Ai = a + r, 

|0> 1 = -^(|0,0> + fl | 
|l>i = ^(|0,l> + e 1 



1.0) ) |2>o = ^(|0,0)-|l,0» (104) 

1.1) ) |3) = ^(|0,1)-|1,1», (105) 

0,1)) 12)00 = ^(|0,0)-|0,1» (106) 

1,1)) = ^(|l,0)+i|l,l». (107) 

then generates the third basis 

1,1)) |2>i = ^(|0,0> + ^ 2 |1, 1» (108) 

1,0)) |3) 1 = ^(|O,1) + 3 |1,O)), (109) 



unbiased with respect to eqs. (105) and (107). We have removed an overall 
phase from each vector leaving the remaining free phases 9q, . . . ,63. The 
conditions for the vectors to form a basis are simply that 1 + 9 9 2 = and 
1 + 9 1 T 3 = 0. 

In dimension N = A, this method constructs the complete set of solutions 
to the cyclic iV-roots problem. In dimension 9 however, whilst we find the 
two parameter family of solutions, there are 6,156 other isolated points [30]. 
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9. Conclusions and questions 

Our main result is that the entire Clifford group admits a representation using 
only monomial phase-permutation matrices if and only if the dimension is a 
square number. We also gave such a representation explicitly. 

In the course of the proof we proved some theorems about the Weyl- 
Heisenberg and Clifford groups that we suspect are known, but which we 
were unable to find in the literature. 

We used this representation to gain some insight into the SIC problem. 
It shares with the standard representation the elegant property that a SIC 
when orthogonally projected onto the simplex spanned by the basis within 
the body of density matrices is again a regular simplex, with iV rather than 
iV 2 vertices. Expressions for the SIC fiducials are much simplified, but are 
still considerably hard to calculate when the dimension exceeds 2 2 . At least 
the case N = 4 is now trivial, and we can now solve the case N = 9 by 
hand. More significantly, we find for the first time an exact solution to the 
SIC problem in dimension N = 16. 

We have applied the phase-permutation representation to the problem of 
finding all vectors unbiased with respect to both the standard and the Fourier 
bases. Families of solutions can then be constructed naturally from Latin 
squares. Finally, we note that quantum mechanics in square-dimensional 
Hilbert spaces is of particular importance because they admit bipartite en- 
tanglement; we therefore expect that the phase-permutation representation 
will have many other applications. 
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Appendix A: Action of the Clifford group 

Lemma 3: The action of the Clifford group on H(N)/Z(N) ~ Z 2 N is iso- 
morphic to SL(2,N). 

Proof. Let G : 1? N — > I? N be a transformation induced by the action of a 
Clifford unitary on H(N)/Z(N). First, we show that G must be an element 
of SL(2,N). 
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That G is linear follows from the fact that H(N)/Z(N) is a projective 
representation of 1? N . Now consider the following commutation relation, 
which is a simple consequence of (8): 



DijDki — uj 



kj - U D kl D tJ 



(110) 



for k, I G Zjv- Conjugate every matrix appearing in the relation above by 
the Clifford unitary Uq. All the phase factors r k appearing in the definition 
(10) cancel, because they occur on both sides of the equality. With (i',f) = 
G(i,j) and (k',l') = G(k,l), we conclude that 



Because co has order N, G preserves symplectic inner products modulo N. 
Thus, G G SL(2, N) as claimed. 

Next, we have to show that every transformation in SL(2,N) can be 
realized. For iV odd, this is the content of (14) proven in Ref. [14]. Hence, 
we only need to consider the case of even N. In this case, (14) says that 
if G G SL(2,N), then GmodiV may be realized as a transformation of 
H(N)/Z(N). Therefore, what remains to be shown is that every matrix G 
in SL(2, N) can be written as G modiV for some G G SL(2, 2N). 

Write N as N = 2 l n for n odd. In Appendix B, we show that 



H(N) ~ H{2 1 ) x H(n), SL(2, N) ~ SL(2, 22 l ) x SL(2, n). (112) 



What is more, "the even and the odd parts do not mix" in the sense that 
SL(2, 22 l ) only acts on H(2 l ) and SL(2, n) only acts on H(n). Therefore, we 
need to prove the claim only for the case N = 2 l . 

So let G G SL(2, N). Then det G = kN + 1 for some integer k. If k is 
even, then det G = 1 mod 2N and therefore G G SL(2, 2N), so we are done. 
Thus we assume that k is odd. Not all matrix elements of G are even, for 
then the range of G would consist only of vectors with even components. 
This would contradict that fact that G is invertible. Assume for now that a, 
the top left matrix element of G, is odd (we label the matrix elements of G 
as in (12)). Then it has an inverse a' 1 modulo 2N . Now let 



UJ 



kj-il _ . ,k'j'-i'l' 



(111) 




(113) 
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Then 

detG = detG + aa- 1 N=(k + aa- 1 )N+l = l mod2iV. (114) 

Thus G E SL(2,2N). The cases where one of the other matrix elements of 
G is odd are treated analogously. □ 

Appendix B: Tensor Product Representation 

The Weyl-Heisenberg group in dimension N is defined in terms of three 
generators X, Z and f and the relations between them. We denote the 
abstract group elements with a tilde; their standard unitary representations, 
defined in eq. (3), appear without. In odd dimensions, the Weyl-Heisenberg 
group, H(N) is given by 

(x,Z,f : X N = Z N = f N = 1, ZX = fXZ, fX = Xf, fZ = Zrj ,(115) 

whilst for N even, we choose to enlarge the centre and define H(N) to be 

^X, Z,f : X N = Z N = f 2N — 1, ZX — t 2 XZ, fX = Xf, fZ = ZfJ . 

(116) 

When the dimension has the prime factorization N = n\ni . . . n r , where 
rij = p U j 3 , the group is a direct product of smaller groups, 

H(N) = H( ni ) x H(n 2 ) x • • • x H(n r ). (117) 

To see this, we construct an isomorphism as follows. Let the group elements 
of H(rij) be generated by Xj, Zj and fj and define elements of H(ni) x • • • x 
H{n r ) as 

x = (Xi, . . . , X r ), 

z = (Z u ...,Z r ), (118) 
t = (fi,...,f r ). 

The elements x, z and t satisfy the relations for the group H(N) since, for 
example, 

X N = ((X^) N/ni , (X™ 2 ) N/n \ (X™ r ) N/nr ) = (1, 1, . . . , 1). (119) 
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Therefore, the map 

9 : X a Z b r c ->■ x a z b t c (120) 

is a homomorphism. 

The image of 6 is given by all elements of the form x a z b t c and we now 
show that it contains the group H(rii) x • • • x H(n r ). The Chinese remainder 
theorem tells us that since rij and n& are coprime for all j ^ k there exists 
an integer Ai such that Ai = 1 mod n\ and Ai = mod rij for j = 2, . . . , r. 
The integer Ai picks out the first component of x, 

x^ = (X?\X*\ X r Al ) = (X u 1, . . . , 1). (121) 

In the same way, there exist integers, A 2 , . . . , A r and /ii, . . . , fj, r such that 

x x > = (1,...,1,1,,1,...,1) 

= (1,...,1,%1,...,1). (122) 

The components of the element t are computed modulo rij so we need 
to modify our argument slightly. In even dimensions, the Chinese remainder 
theorem still applies since only one of the factors, say rii, is even. The 
integer 2n\ is therefore coprime to rij for all j, and we can again find integers 
vi, . . . , v r such that 

I" (1 1!. (123) 

Finally, the size of the two groups are equal, \H(N)\ = N 2 N = \H(ni) x 
• • • x H(n r )\ so 6 is an isomorphism. 

Now for the Clifford group C(N). We use the fact that C(N) is the 
semi-direct product of H(N) and SL(2,N). We start by taking elements of 
SL(2,N) and computing their components modulo flj, that is, 




(124) 



where aij = a mod rij, f3j = (5 mod rij, 7j = 7 mod rij and 5j = 5 mod rij. 
Then the map 

T : SL(2, N) ->■ SL(2, n x ) x • • • x 5L(2, n r ), (125) 

defined by 

T(F) = (F 1 ,...,F r ), (126) 
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is an isomorphism. The proof follows a similar argument to the above and 
implies that 

C(N) ~ H(m) x ••• x H(n r ) x SL(2,fn) x ••• x SL(2,n r ) 

~ if (rii) x SX(2, ni) x • • • x #(n r ) x SX(2, n r ) (127) 
~ C(ni) x • • • x C(n r ). 

These two observations mean that every displacement operator can be 
written as a tensor product of displacement operators in smaller Hilbert 
spaces because of the following fact from finite group theory (see for example 
Theorem 19.18 of Ref. [31]). Let G and J be groups, then every irreducible 
representation of the group G x J is a tensor product of an irreducible repre- 
sentation of G with an irreducible representation of J. The Weyl-Heisenberg 
and Clifford groups are direct products of the groups defined over the prime 
factorization and therefore all irreducible representations can be written as 
a tensor product of irreducible representations of the smaller groups. 

Care is required when writing down the isomorphisms in terms of the 
standard unitary representation, X, Z and r, defined in eq. (3). If we define 
the map 

r] : x a z b t c ->• (rf . . . r r c ) {X\Z\ <g> • • • <g> X^Z^) , (128) 
we have the problem that in general, 

Tf.-.T^T 6 , (129) 

meaning that the right hand side of eq. (128) cannot be the image of 9 under 
any unitary induced mapping. To fix this problem, we redefine i] to be 

rj' : x a z b t c -)• (rf lC . . . r^ c ) (X?Z« lb <g> • • • <g> X a r Z^ h ) , (130) 

where Kj is the multiplicative inverse of N/rij mod fij. The map r( then 
satisfies all of the required properties to be an isomorphism. 

To construct the isomorphism for the standard unitary representation of 
the Clifford group, we take the image of the symplectic matrices, Fj to be 
Uf', where 

3 

f; = ( Z> t ft ) ■ (131) 

rather than Upy Whilst we did not need the explicit form of the two isomor- 
phisms in this paper, we hope that it will prove a useful tool elsewhere. 



40 



Appendix C: Fiducial vectors for N = 16 



16a 

w2:=Sqrt(2) ; 
wl3:=-Sqrt(13) ; 
wl7:=-Sqrt(17) ; 

I:=Sqrt(-l); 

w26:=w2*wl3; 

w34:=w2*wl7; 
w221:=wl3*wl7; 
w442:=w2*wl3*wl7; 

r2:=Sqrt(wl3*wl7-ll) ; 
r3:=Sqrt(15+wl7) ; 

tl:=Sqrt(15+(4-wl7)*r3-3*wl7) ; 

t2 : =Sqrt ( ( ( (-5*wl7+3) *wl3+ (39*wl7-65) ) *r3+ ( ( 16*wl7-72) *wl3+936) ) *t l-208*wl3+2288) ; 
t3:=Sqrt(2-w2) ; 
t4:=Sqrt(2+t3) ; 



c01:=(-20/13*w26*r3*tl*t2+((-90*w2+60*wl3-20*wl7+30*w26+10*w34+20*w221+10*w442-60)*r3 

+(60*wl3+60*wl7+20*w221-140))*tl*t3+((350*w2+160*wl3-40*wl7+150*w26-50*w34-10*w442+760)*r3 

+(400*w2+440*wl3-600*wl7+320*w26-320*w34-40*w221-80*w442+360))*t3)*I+(-40/13*wl3 

-20/13*w26)*r3*tl*t2+((140*w2+20*w26+160)*r3-20*w2-260*wl3-500*wl7-140*w26-140*w34-60*w221 

-20*w442-1260)*tl*t3+((530*w2-60*wl3-100*wl7-30*w26-70*w34+20*w221+10*w442+620)*r3+(520*w2 

-400*wl3-1120*wl7-200*w26-520*w34+40*w442-1200))*t3; 

c02 : =( ( (10*w2-4*w26-2*w34) *r2*r3+(-24*w2-24*w26-12*w34-4*w442) *r2+(-50*wl3+10*wl7-10*w221 

-30)*r3-20*wl3-140*wl7-20*w221-540)*tl*t3*t4+((-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221 

-3/2*w442-ll)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3 

+5*»17+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((-10*w2-27*wl3-21*wl7+4*w26+2*w34-3*w221-53)*r2*r3+(24*w2-24*wl3 

+8*wl7+24*w26+12*w34-8*w221+4*w442-88)*r2+(-85*w2+50*wl3-10*wl7-25*w26+15*w34+10*w221 

-5*w442+30)*r3+(80*w2+20*wl3+140*wl7+40*w26+40*w34+20*w221+540))*tl*t4+((53/2*w2+ll*wl3 

+3*wl7+27/2*w26+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34 

+14*w221+4*w442+174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((-36*w2-28*w26-4*w34+4*w442)*r2*r3+(122*w2-86*w26+22*w34 

+22*w442)*r2+(-80*wl3+40*wl7-600)*r3-220*wl3+260*wl7-20*w221-1300)*t3*t4+((-46*w2+25*wl3 

+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34-18*w221 

-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((36*w2-56*wl3-28*wl7+28*w26+4*w34-4*w221-4*w442 

-144)*r2*r3+(-122*w2-64*wl3+48*wl7+86*w26-22*w34+16*w221-22*w442+96)*r2+(-470*w2+80*wl3 

-40*wl7+10*w26+30*w34-10*w442+600)*r3-240*w2+220*wl3-260*wl7+40*w26+400*w34+20*w221 

-40*w442+1300) *t4+( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3+(-48*w2 

-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (-3*wl3 

-9*wl7+w221+ll)*r2*r3+(56*wl3+48*wl7+24*w221+304)*r2+(-55*w2+5*w26-15*w34+5*w442)*r3 

-60*w2+100*w26+100*w34+20*w442)*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7+37/2*w26+31/2*w34+w221 

+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221-3*w442-132)*r2+(55*w2 

-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((-32*w2+3*wl3+9*wl7-14*w26-12*w34-w221-2*w442-ll)*r2*r3 



41 



+(130*w2-56*wl3-48*wl7+42*w26+26*w34-24*w221+10*w442-304)*r2+(55*w2+10*wl3-10*wl7-5*w26 

+15*w34+10*w221-5*w442-90)*r3+(60*w2+40*wl3-80*wl7-100*w26-100*w34-20*w442-680))*tl*t4 

+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2-42*wl3 

-26*wl7+12*w26+16*w34-10*w221-130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442 

+55)*r3-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620)*tl+((-46*wl3-38*wl7+6*w221 

+86) *r2*r3+ (56*wl3+8*wl7+72*w221+952) *r2+ (-60*w2+20*w26+40*w34) *r3-280*w2+160*w26 

+120*w34) *t3*t4+( (75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442-7) *r2*r3+ (-130*w2 

-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34 

+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3 

+((-42*w2+46*wl3+38*wl7-26*w26-18*w34-6*w221-2*w442-86)*r2*r3+(286*w2-56*wl3-8*wl7+62*w26 

+26*w34-72*w221+26*w442-952) *r2+ (60*w2+80*wl3-20*w26-40*w34+160) *r3+ (280*w2+220*wl3 

-340*wl7-160*w26-120*w34-60*w221-940))*t4+(-36*w2+26*wl3+18*wl7-10*w34+2*w221-6*w442 

+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221-14*w442-286)*r2+(80*w2-90*wl3 

+10*wl7+40*w26+10*w221-130)*r3-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442 

-1060; 

c03:=(((-ll/26*wl3-l/2*wl7-3/26*w221-l/2)*r2*r3+(10*w2+20/13*w26+10/13*w442)*r2 

-40/13*wl3*r3)*tl*t2+((-43*w2-27*wl3-21*wl7-ll*w26-3*w34-3*w221-3*w442-53)*r2*r3+(-174*w2 

-24*wl3+8*wl7-14*w26-22*w34-8*w221-14*w442-88)*r2+(-55*w2-10*wl3+10*wl7+25*w26+5*w34 

-10*w221+5*w442+90) *r3+ (620*w2-40*wl3+80*wl7+60*w26+180*w34+20*w442+680) ) *t 1+ ( (-128*w2 

-56*wl3-28*wl7+20*w26+20*w34-4*w221-8*w442-144)*r2*r3+(-666*w2-64*wl3+48*wl7+6*w26+18*w34 

+16*w221-46*w442+96)*r2+(130*w2-80*wl3+90*w26-10*w34-10*w442-160)*r3+(1060*w2-220*wl3 

+340*wl7+260*w26+140*w34+60*w221-20*w442+940) ) ) *I+ ( (H/26*wl3+l/2*wl7+3/26*w221+l/2) *r2*r3 

+(10*w2+20/13*w26+10/13*w442)*r2)*tl*t2+((-32*w2+19*wl3+17*wl7-14*w26-12*w34+3*w221-2*w442 

+73)*r2*r3+(130*w2-24*wl3-32*wl7+42*w26+26*w34+10*w442+40)*r2+(-85*w2+20*wl7-25*w26 

+15*w34-5*w442+20)*r3+(80*w2-160*wl3-280*wl7+40*w26+40*w34-40*w221-560))*tl+(-42*w2 

+20*wl7-26*w26-18*w34+12*w221-2*w442+72)*r2*r3+(286*w2-108*wl3-4*wl7+62*w26+26*w34 

+28*w221+26*w442+148) *r2+ (-470*w2+40*wl3-80*wl7+10*w26+30*w34-10*w442+280) *r3-240*w2 

-100*wl3-580*wl7+40*w26+400*w34-60*w221-40*w442-380; 

c04 : = ( ( (-63*w2-23*w26-19*w34-3*w442) *r2*r3+ (-64*w2+20*w34-4*w442) *r2+ (-20*wl7+140) *r3 

-100*wl3-220*wl7-20*w221-700)*tl*t3*t4+((53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221+3/2*w442 

+Il)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3+5*wl7 

+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221-10*w442 

-10)*tl*t3+((63*w2-19*wl3-17*wl7+23*w26+19*w34-3*w221+3*w442-73)*r2*r3+(64*w2+24*wl3 

+32*wl7-20*w34+4*w442-40)*r2+(55*w2+20*wl7-25*w26-5*w34-5*w442-140)*r3-620*w2+100*wl3 

+220*wl7-60*w26-180*w34+20*w221-20*w442+700)*tl*t4+((53/2*w2+ll*wl3+3*wl7+27/2*w26 

+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34+14*w221+4*w442 

+174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26+140*w34 

+20*w442-80) ) *t 1+ ( (-108*w2-28*w26-24*w34-8*w442) *r2*r3+ (-26*w2+22*w26+26*w34-6*w442) *r2 

+ (-100*wl3-20*wl7+20*w221+340) *r3-300*wl3-540*wl7+60*w221-820) *t3*t4+ ( (46*w2-25*wl3 

-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 

+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((108*w2-20*wl7+28*w26+24*w34-12*w221+8*w442 

-72)*r2*r3+(26*w2+108*wl3+4*wl7-22*w26-26*w34-28*w221+6*w442-148)*r2+(-130*w2+100*wl3 

+20*wl7-90*w26+10*w34-20*w221+10*w442-340)*r3-1060*w2+300*wl3+540*wl7-260*w26-140*w34 

-60*w221+20*w442+820) *t4+ ( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3 

+(-48*w2-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(-140*w2-10*wl3-30*wl7-20*w26 

+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221+30*w442+240)))*I 

+ ( (-25*wl3-15*wl7-5*w221-75) *r2*r3+ (28*wl3+4*wl7-4*w221-44) *r2+ (35*w2-5*w26-5*w34 

-5*w442)*r3+(620*w2+60*w26+180*w34+20*w442))*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26 

-31/2*w34-w221-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442 

+132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7 

-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+((-43*w2+25*wl3+15*wl7-ll*w26-3*w34+5*w221 

-3*w442+75)*r2*r3+(-174*w2-28*wl3-4*wl7-14*w26-22*w34+4*w221-14*w442+44)*r2+(-35*w2 
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-20*wl7+5*w26+5*w34+5*w442-20)*r3-620*w2+160*wl3+280*wl7-60*w26-180*w34+40*w221-20*w442 

+560)*tl*t4+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2 

-42*wl3-26*wl7+12*w26+16*w34-10*w221-130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 

-5*w442-55)*r3+(340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620))*tl+((-6*wl3+2*wl7 

-10*w221-170)*r2*r3+(68*wl3+44*wl7-20*w221-380)*r2+(-220*w2-60*w26+40*w34)*r3+(660*w2 

-60*w26+460*w34+60*w442))*t3*t4+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442 

+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3 

-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34 

+60*w221+10*w442+60)*t3+((-128*w2+6*wl3-2*wl7+20*w26+20*w34+10*w221-8*w442+170)*r2*r3 

+ (-666*w2-68*wl3-44*wl7+6*w26+18*w34+20*w221-46*w442+380) *r2+ (220*w2-40*wl3+80*wl7+60*w26 

-40*w34-280) *r3-660*w2+100*wl3+580*wl7+60*w26-460*w34+60*w221-60*w442+380) *t4+ (-36*w2 

+26*wl3+18*wl7-10*w34+2*w221-6*w442+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221 

-14*w442-286)*r2+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7 

-110*w26+170*w34-20*w221+30*w442+1060; 

c05:=(-20/13*w26*r3*tl*t2+((20*w2+30*wl3+30*wl7+12*w26+16*w34+2*w221+42)*r2*r3+(-110*w2 

-32*wl3-56*wl7-14*w26-42*w34-16*w221-10*w442-216)*r2)*tl*t3+((-20*w2+102*wl3+66*wl7+48*w26 

+44*w34-2*w221+58)*r2*r3+(-640*w2+8*wl3-56*wl7-16*w26-8*w34-88*w221-40*w442 

-1048) *r2) *t3) *I+(40/13*wl3-20/13*w26) *r3*tl*t2+( (42*w2+44*wl3+32*wl7+10*w26+10*w34+8*w221 

+2*w442+148)*r2*r3+(-154*w2-52*wl3-36*wl7-66*w26-38*w34+4*w221-14*w442+84)*r2)*tl*t3 

+ ( (6*w2+6*wl3+18*wl7-2*w26+14*w34+22*w221+6*w442+242) *r2*r3+ (-164*w2-176*wl3-48*wl7 

-148*w26-4*w34+48*w221-4*w442+528)*r2)*t3; 

c06 : = ( ( (-63*w2-23*w26-19*w34-3*w442) *r2*r3+ (-64*w2+20*w34-4*w442) *r2+ (20*wl7-140) *r3 

+(100*wl3+220*wl7+20*w221+700))*tl*t3*t4+((-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221 

-3/2*w442-ll)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*»17-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((63*w2-19*wl3-17*wl7+23*w26+19*w34-3*w221+3*w442-73)*r2*r3+(64*w2 

+24*wl3+32*wl7-20*w34+4*w442-40)*r2+(-55*w2-20*wl7+25*w26+5*w34+5*w442+140)*r3+(620*w2 

-100*wl3-220*wl7+60*w26+180*w34-20*w221+20*w442-700))*tl*t4+((-53/2*w2-ll*wl3-3*wl7 

-27/2*w26-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34-14*w221 

-4*w442-174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26 

+140*w34+20*w442-80) ) *t 1+ ( (-108*w2-28*w26-24*w34-8*w442) *r2*r3+ (-26*w2+22*w26+26*w34 

-6*w442) *r2+ (100*wl3+20*wl7-20*w221-340) *r3+ (300*wl3+540*wl7-60*w221+820) ) *t3*t4+ ( (-46*w2 

+25*wl3+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34 

-18*w221-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3 

-420*wl7+10*w26-430*w34-20*w221-10*w442+460)*t3+((108*w2-20*wl7+28*w26+24*w34-12*w221 

+8*w442-72)*r2*r3+(26*w2+108*wl3+4*wl7-22*w26-26*w34-28*w221+6*w442-148)*r2+(130*w2 

-100*wl3-20*wl7+90*w26-10*w34+20*w221-10*w442+340)*r3+(1060*w2-300*wl3-540*wl7+260*w26 

+140*w34+60*w221-20*w442-820) ) *t4+( (-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442 

-128)*r2*r3+(48*w2+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(-140*w2-10*wl3 

-30*wl7-20*w26+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221 

+30*w442+240) ) ) *I+ ( (-25*wl3-15*wl7-5*w221-75) *r2*r3+ (28*wl3+4*wl7-4*w221-44) *r2+ (-35*w2 

+5*w26+5*w34+5*w442)*r3-620*w2-60*w26-180*w34-20*w442)*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7 

+37/2*w26+31/2*w34+w221+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221 

-3*w442-132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3 

-70*wl7-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+((-43*w2+25*wl3+15*wl7-ll*w26-3*w34 

+5*w221-3*w442+75)*r2*r3+(-174*w2-28*wl3-4*wl7-14*w26-22*w34+4*w221-14*w442+44)*r2+(35*w2 

+20*wl7-5*w26-5*w34-5*w442+20)*r3+(620*w2-160*wl3-280*wl7+60*w26+180*w34-40*w221+20*w442 

-560))*tl*t4+((73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32)*r2*r3+(20*w2 

+42*wl3+26*wl7-12*w26-16*w34+10*w221+130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 

-5*w442-55)*r3+(340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620))*tl+((-6*wl3+2*»17 

-10*w221-170)*r2*r3+(68*wl3+44*wl7-20*w221-380)*r2+(220*w2+60*w26-40*w34)*r3-660*w2+60*w26 

-460*w34-60*w442) *t3*t4+ ( (75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442-7) *r2*r3 

+ (-130*w2-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402) *r2+ (-35*w2-90*wl3-10*wl7-55*w26 
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-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442 

+60)*t3+((-128*w2+6*wl3-2*wl7+20*w26+20*w34+10*w221-8*w442+170)*r2*r3+(-666*w2-68*wl3 

-44*wl7+6*w26+18*w34+20*w221-46*w442+380)*r2+(-220*w2+40*wl3-80*wl7-60*w26+40*w34+280)*r3 

+(660*w2-100*wl3-580*wl7-60*w26+460*w34-60*w221+60*w442-380))*t4+(36*w2-26*wl3-18*wl7 

+10*w34-2*w221+6*w442-42) *r2*r3+ (74*w2+62*wl3+26*wl7-54*w26-2*w34+26*w221+14*w442+286) *r2 

+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7-110*w26+170*w34 

-20*w221+30*w442+1060 ; 

c07 :=( ( (H/26*wl3+l/2*wl7+3/26*w221+l/2) *r2*r3+ (-10*w2-20/13*w26-10/13*w442) *r2 

-40/13*wl3*r3)*tl*t2+((-43*w2-27*wl3-21*wl7-ll*w26-3*w34-3*w221-3*w442-53)*r2*r3+(-174*w2 

-24*wl3+8*wl7-14*w26-22*w34-8*w221-14*w442-88)*r2+(55*w2+10*wl3-10*wl7-25*w26-5*w34 

+10*w221-5*w442-90) *r3-620*w2+40*wl3-80*wl7-60*w26-180*w34-20*w442-680) *tl+ ( (-128*w2 

-56*wl3-28*wl7+20*w26+20*w34-4*w221-8*w442-144)*r2*r3+(-666*w2-64*wl3+48*wl7+6*w26+18*w34 

+16*w221-46*w442+96)*r2+(-130*w2+80*wl3-90*w26+10*w34+10*w442+160)*r3-1060*w2+220*wl3 

-340*wl7-260*w26-140*w34-60*w221+20*w442-940))*I+((-ll/26*wl3-l/2*wl7-3/26*w221-l/2)*r2*r3 

+(-10*w2-20/13*w26-10/13*w442)*r2)*tl*t2+((-32*w2+19*wl3+17*wl7-14*w26-12*w34+3*w221-2*w442 

+73)*r2*r3+(130*w2-24*wl3-32*wl7+42*w26+26*w34+10*w442+40)*r2+(85*w2-20*wl7+25*w26-15*w34 

+5*w442-20)*r3-80*w2+160*wl3+280*wl7-40*w26-40*w34+40*w221+560)*tl+(-42*w2+20*wl7-26*w26 

-18*w34+12*w221-2*w442+72)*r2*r3+(286*w2-108*wl3-4*wl7+62*w26+26*w34+28*w221+26*w442 

+148) *r2+ (470*w2-40*wl3+80*wl7-10*w26-30*w34+10*w442-280) *r3+240*w2+100*wl3+580*wl7 

-40*w26-400*w34+60*w221+40*w442+380; 

c08:=(((10*w2-4*w26-2*w34)*r2*r3+(-24*w2-24*w26-12*w34-4*w442)*r2+(50*wl3-10*wl7+10*w221 

+30) *r3+ (20*wl3+140*wl7+20*w221+540) ) *t I*t3*t4+ ( (53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221 

+3/2*w442+ll)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((-10*w2-27*wl3-21*wl7+4*w26+2*w34-3*w221-53)*r2*r3+(24*w2-24*wl3 

+8*wl7+24*w26+12*w34-8*w221+4*w442-88)*r2+(85*w2-50*wl3+10*wl7+25*w26-15*w34-10*w221 

+5*w442-30)*r3-80*w2-20*wl3-140*wl7-40*w26-40*w34-20*w221-540)*tl*t4+((-53/2*w2-ll*wl3 

-3*wl7-27/2*w26-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34 

-14*w221-4*w442-174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((-36*w2-28*w26-4*w34+4*w442)*r2*r3+(122*w2-86*w26+22*w34 

+22*w442)*r2+(80*wl3-40*wl7+600)*r3+(220*wl3-260*wl7+20*w221+1300))*t3*t4+((46*w2-25*wl3 

-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 

+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((36*w2-56*wl3-28*wl7+28*w26+4*w34-4*w221-4*w442 

-144) *r2*r3+ (-122*w2-64*wl3+48*wl7+86*w26-22*w34+16*w221-22*w442+96) *r2+ (470*w2-80*wl3 

+40*wl7-10*w26-30*w34+10*w442-600)*r3+(240*w2-220*wl3+260*wl7-40*w26-400*w34-20*w221 

+40*w442-1300))*t4+((-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442-128)*r2*r3+(48*w2 

+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (-3*wl3 

-9*wl7+w221+ll)*r2*r3+(56*wl3+48*wl7+24*w221+304)*r2+(55*w2-5*w26+15*w34-5*w442)*r3 

+(60*w2-100*w26-100*w34-20*w442))*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26-31/2*w34-w221 

-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442+132)*r2+(55*w2 

-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((-32*w2+3*wl3+9*wl7-14*w26-12*w34-w221-2*w442-ll)*r2*r3 

+(130*w2-56*wl3-48*wl7+42*w26+26*w34-24*w221+10*w442-304)*r2+(-55*w2-10*wl3+10*wl7+5*w26 

-15*w34-10*w221+5*w442+90)*r3-60*w2-40*wl3+80*wl7+100*w26+100*w34+20*w442+680)*tl*t4 

+ ( (73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32) *r2*r3+ (20*w2+42*wl3+26*wl7 

-12*w26-16*w34+10*w221+130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442+55)*r3 

-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620) *tl+( (-46*wl3-38*wl7+6*w221+86) *r2*r3 

+(56*wl3+8*wl7+72*w221+952) *r2+(60*w2-20*w26-40*w34) *r3+(280*w2-160*w26-120*w34) ) *t3*t4 

+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26 

+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3 

-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3+((-42*w2+46*wl3+38*wl7 
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-26*w26-18*w34-6*w221-2*w442-86)*r2*r3+(286*w2-56*wl3-8*wl7+62*w26+26*w34-72*w221+26*w442 

-952)*r2+(-60*w2-80*wl3+20*w26+40*w34-160)*r3-280*w2-220*wl3+340*wl7+160*w26+120*w34 

+60*w221+940)*t4+(36*w2-26*wl3-18*wl7+10*w34-2*w221+6*w442-42)*r2*r3+(74*w2+62*wl3 

+26*wl7-54*w26-2*w34+26*w221+14*w442+286)*r2+(80*w2-90*wl3+10*wl7+40*w26+10*w221-130)*r3 

-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442-1060; 

c09:=(-20/13*w26*r3*tl*t2+((90*w2-60*wl3+20*wl7-30*w26-10*w34-20*w221-10*w442+60)*r3-60*wl3 

-60*wl7-20*w221+140)*tl*t3+((-350*w2-160*wl3+40*wl7-150*w26+50*w34+10*w442-760)*r3-400*w2 
-440*wl3+600*wl7-320*w26+320*w34+40*w221+80*w442-360)*t3)*I+(-40/13*wl3-20/13*w26)*r3*tl*t2 
+ ( (-140*w2-20*w26-160) *r3+ (20*w2+260*wl3+500*wl7+140*w26+140*w34+60*w221+20*w442 
+1260))*tl*t3+((-530*w2+60*wl3+100*wl7+30*w26+70*w34-20*w221-10*w442-620)*r3-520*w2 
+400*wl3+1120*wl7+200*w26+520*w34-40*w442+1200)*t3; 

clO : = ( ( (-10*w2+4*w26+2*w34) *r2*r3+ (24*w2+24*w26+12*w34+4*w442) *r2+ (50*wl3-10*wl7+10*w221 

+30) *r3+ (20*wl3+140*wl7+20*w221+540) ) *t I*t3*t4+ ( (-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221 

-3/2*w442-ll)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((10*w2+27*wl3+21*wl7-4*w26-2*w34+3*w221+53)*r2*r3+(-24*w2+24*wl3 

-8*wl7-24*w26-12*w34+8*w221-4*w442+88)*r2+(85*w2-50*wl3+10*wl7+25*w26-15*w34-10*w221 

+5*w442-30)*r3-80*w2-20*wl3-140*wl7-40*w26-40*w34-20*w221-540)*tl*t4+((53/2*w2+ll*wl3 

+3*ul7+27/2*w26+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34 

+14*w221+4*w442+174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((36*w2+28*w26+4*w34-4*w442)*r2*r3+(-122*w2+86*w26-22*w34 

-22*w442) *r2+ (80*wl3-40*wl7+600) *r3+ (220*wl3-260*wl7+20*w221+1300) ) *t3*t4+ ( (-46*w2+25*wl3 

+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34-18*w221 

-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((-36*w2+56*wl3+28*wl7-28*w26-4*w34+4*w221+4*w442 

+144)*r2*r3+(122*w2+64*wl3-48*wl7-86*w26+22*w34-16*w221+22*w442-96)*r2+(470*w2-80*wl3 

+40*wl7-10*w26-30*w34+10*w442-600)*r3+(240*w2-220*wl3+260*wl7-40*w26-400*w34-20*w221 

+40*w442-1300) ) *t4+ ( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3+ (-48*w2 

-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (3*wl3 

+9*wl7-w221-ll)*r2*r3+(-56*wl3-48*wl7-24*w221-304)*r2+(55*w2-5*w26+15*w34-5*w442)*r3 

+(60*w2-100*w26-100*w34-20*w442))*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7+37/2*w26+31/2*w34+w221 

+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221-3*w442-132)*r2+(55*w2 

-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((32*w2-3*wl3-9*wl7+14*w26+12*w34+w221+2*w442+ll)*r2*r3 

+(-130*w2+56*wl3+48*wl7-42*w26-26*w34+24*w221-10*w442+304)*r2+(-55*w2-10*wl3+10*wl7+5*w26 

-15*w34-10*w221+5*w442+90)*r3-60*w2-40*wl3+80*wl7+100*w26+100*w34+20*w442+680)*tl*t4 

+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2-42*wl3 

-26*wl7+12*w26+16*w34-10*w221-130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442 

+55)*r3-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620)*tl+((46*wl3+38*wl7-6*w221 

-86)*r2*r3+(-56*wl3-8*wl7-72*w221-952)*r2+(60*w2-20*w26-40*w34)*r3+(280*w2-160*w26 

-120*w34))*t3*t4+((75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442-7)*r2*r3+(-130*w2 

-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34 

+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3 

+ ( (42*w2-46*wl3-38*wl7+26*w26+18*w34+6*w221+2*w442+86) *r2*r3+ (-286*w2+56*wl3+8*wl7-62*w26 

-26*w34+72*w221-26*w442+952)*r2+(-60*w2-80*wl3+20*w26+40*w34-160)*r3-280*w2-220*wl3 

+340*wl7+160*w26+120*w34+60*w221+940)*t4+(-36*w2+26*wl3+18*wl7-10*w34+2*w221-6*w442 

+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221-14*w442-286)*r2+(80*w2-90*wl3 

+ 10*wl7+40*w26+10*w221-130)*r3-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442 

-1060; 

cll:=(((-ll/26*wl3-l/2*wl7-3/26*w221-l/2)*r2*r3+(10*w2+20/13*w26+10/13*w442)*r2 

-40/13*wl3*r3)*tl*t2+((43*w2+27*wl3+21*wl7+ll*w26+3*w34+3*w221+3*w442+53)*r2*r3+(174*w2 



45 



+24*wl3-8*wl7+14*w26+22*w34+8*w221+14*w442+88)*r2+(55*w2+10*wl3-10*wl7-25*w26-5*w34 

+10*w221-5*w442-90) *r3-620*w2+40*wl3-80*wl7-60*w26-180*w34-20*w442-680) *t 1+ ( ( 128*w2 

+56*wl3+28*wl7-20*w26-20*w34+4*w221+8*w442+144)*r2*r3+(666*w2+64*wl3-48*wl7-6*w26-18*w34 

-16*w221+46*w442-96)*r2+(-130*w2+80*wl3-90*w26+10*w34+10*w442+160)*r3-1060*w2+220*wl3 

-340*wl7-260*w26-140*w34-60*w221+20*w442-940))*I+((ll/26*wl3+l/2*wl7+3/26*w221+l/2)*r2*r3 

+(10*w2+20/13*w26+10/13*w442)*r2)*tl*t2+((32*w2-19*wl3-17*wl7+14*w26+12*w34-3*w221+2*w442 

-73)*r2*r3+(-130*w2+24*wl3+32*wl7-42*w26-26*w34-10*w442-40)*r2+(85*w2-20*wl7+25*w26 

-15*w34+5*w442-20)*r3-80*w2+160*wl3+280*wl7-40*w26-40*w34+40*w221+560)*tl+(42*w2-20*wl7 

+26*w26+18*w34-12*w221+2*w442-72)*r2*r3+(-286*w2+108*wl3+4*wl7-62*w26-26*w34-28*w221 

-26*w442-148)*r2+(470*w2-40*wl3+80*wl7-10*w26-30*w34+10*w442-280)*r3+240*w2+100*wl3 

+580*wl7-40*w26-400*w34+60*w221+40*w442+380; 

Cl2:=(((63*w2+23*w26+19*w34+3*w442)*r2*r3+(64*w2-20*w34+4*w442)*r2+(20*wl7-140)*r3+(100*wl3 
+220*wl7+20*w221+700))*tl*t3*t4+((53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221+3/2*w442 
+Il)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3+5*wl7 
+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221-10*w442 
-10)*tl*t3+((-63*w2+19*wl3+17*wl7-23*w26-19*w34+3*w221-3*w442+73)*r2*r3+(-64*w2-24*wl3 
-32*wl7+20*w34-4*w442+40)*r2+(-55*w2-20*wl7+25*w26+5*w34+5*w442+140)*r3+(620*w2-100*wl3 
-220*wl7+60*w26+180*w34-20*w221+20*w442-700))*tl*t4+((53/2*w2+ll*wl3+3*wl7+27/2*w26 
+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34+14*w221+4*w442 
+174) *r2+ (-10*w2+25*wl3-15*wl7-10*w34+5*w221+85) *r3+ (280*w2-40*wl3-40*wl7+80*w26+140*w34 
+20*w442-80))*tl+((108*w2+28*w26+24*w34+8*w442)*r2*r3+(26*w2-22*w26-26*w34+6*w442)*r2 
+(100*wl3+20*wl7-20*w221-340)*r3+(300*wl3+540*wl7-60*w221+820))*t3*t4+( (46*w2-25*wl3 
-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 
+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 
+10*w26-430*w34-20*w221-10*w442+460)*t3+((-108*w2+20*wl7-28*w26-24*w34+12*w221-8*w442 
+72)*r2*r3+(-26*w2-108*wl3-4*wl7+22*w26+26*w34+28*w221-6*w442+148)*r2+(130*w2-100*wl3 
-20*wl7+90*w26-10*w34+20*w221-10*w442+340)*r3+(1060*w2-300*wl3-540*wl7+260*w26+140*w34 
+60*w221-20*w442-820) ) *t4+ ( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3 
+(-48*w2-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(-140*w2-10*wl3-30*wl7-20*w26 
+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221+30*w442+240)))*I 
+((25*wl3+15*wl7+5*w221+75)*r2*r3+(-28*wl3-4*wl7+4*w221+44)*r2+(-35*w2+5*w26+5*w34 
+5*w442)*r3-620*w2-60*w26-180*w34-20*w442)*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26 
-31/2*w34-w221-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442 
+132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7 
-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+((43*w2-25*wl3-15*wl7+ll*w26+3*w34-5*w221 
+3*w442-75)*r2*r3+(174*w2+28*wl3+4*wl7+14*w26+22*w34-4*w221+14*w442-44)*r2+(35*w2+20*wl7 
-5*w26-5*w34-5*w442+20)*r3+(620*w2-160*wl3-280*wl7+60*w26+180*w34-40*w221+20*w442 
-560))*tl*t4+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2 
-42*wl3-26*wl7+12*w26+16*w34-10*w221-130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 
-5*w442-55) *r3+ (340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620) ) *tl+ ( (6*wl3-2*wl7 
+10*w221+170) *r2*r3+ (-68*wl3-44*wl7+20*w221+380) *r2+ (220*w2+60*w26-40*w34) *r3-660*w2 
+60*w26-460*w34-60*w442)*t3*t4+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442 
+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3 
-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34 
+60*w221+10*w442+60)*t3+((128*w2-6*wl3+2*wl7-20*w26-20*w34-10*w221+8*w442-170)*r2*r3 
+(666*w2+68*wl3+44*wl7-6*w26-18*w34-20*w221+46*w442-380)*r2+(-220*w2+40*wl3-80*wl7-60*w26 
+40*w34+280)*r3+(660*w2-100*wl3-580*wl7-60*w26+460*w34-60*w221+60*w442-380))*t4+(-36*w2 
+26*wl3+18*wl7-10*w34+2*w221-6*w442+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221 
-14*w442-286)*r2+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7 
-110*w26+170*w34-20*w221+30*w442+1060; 

cl3:=(-20/13*w26*r3*tl*t2+((-20*w2-30*wl3-30*wl7-12*w26-16*w34-2*w221-42)*r2*r3+(110*w2 

+32*wl3+56*wl7+14*w26+42*w34+16*w221+10*w442+216)*r2)*tl*t3+((20*w2-102*wl3-66*wl7-48*w26 
-44*w34+2*w221-58)*r2*r3+(640*w2-8*wl3+56*wl7+16*w26+8*w34+88*w221+40*w442 
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+1048) *r2) *t3) *I+(40/13*wl3-20/13*w26) *r3*tl*t2+( (-42*w2-44*wl3-32*wl7-10*w26-10*w34-8*w221 
-2*w442-148)*r2*r3+(154*w2+52*wl3+36*wl7+66*w26+38*w34-4*w221+14*w442-84)*r2)*tl*t3 
+((-6*w2-6*wl3-18*wl7+2*w26-14*w34-22*w221-6*w442-242)*r2*r3+(164*w2+176*wl3+48*wl7 
+ 148*w26+4*w34-48*w221+4*w442-528) *r2) *t3 ; 

cl4 : = ( ( (63*w2+23*w26+19*w34+3*w442) *r2*r3+(64*w2-20*w34+4*w442) *r2+(-20*wl7+140) *r3-100*wl3 
-220*wl7-20*w221-700)*tl*t3*t4+((-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221-3/2*w442 
-H)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3+5*wl7 
+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221-10*w442 
-10)*tl*t3+((-63*w2+19*wl3+17*wl7-23*w26-19*w34+3*w221-3*w442+73)*r2*r3+(-64*w2-24*wl3 
-32*wl7+20*w34-4*w442+40)*r2+(55*w2+20*wl7-25*w26-5*w34-5*w442-140)*r3-620*w2+100*wl3 
+220*wl7-60*w26-180*w34+20*w221-20*w442+700)*tl*t4+((-53/2*w2-ll*wl3-3*wl7-27/2*w26 
-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34-14*w221-4*w442 
-174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26+140*w34 
+20*w442-80))*tl+((108*w2+28*w26+24*w34+8*w442)*r2*r3+(26*w2-22*w26-26*w34+6*w442)*r2 
+(-100*wl3-20*wl7+20*w221+340)*r3-300*wl3-540*wl7+60*w221-820)*t3*t4+((-46*w2+25*wl3 
+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34-18*w221 
-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 
+10*w26-430*w34-20*w221-10*w442+460)*t3+((-108*w2+20*wl7-28*w26-24*w34+12*w221-8*w442 
+72)*r2*r3+(-26*w2-108*wl3-4*wl7+22*w26+26*w34+28*w221-6*w442+148)*r2+(-130*w2+100*wl3 
+20*wl7-90*w26+10*w34-20*w221+10*w442-340)*r3-1060*w2+300*wl3+540*wl7-260*w26-140*w34 
-60*w221+20*w442+820)*t4+((-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442-128)*r2*r3 
+(48*w2+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(-140*w2-10*wl3-30*wl7-20*w26 
+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221+30*w442+240)))*I 
+ ( (25*wl3+15*wl7+5*w221+75) *r2*r3+ (-28*wl3-4*wl7+4*w221+44) *r2+ (35*w2-5*w26-5*w34 
-5*w442)*r3+(620*w2+60*w26+180*w34+20*w442))*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7+37/2*w26 
+31/2*w34+w221+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221-3*w442 
-132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7 
-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+((43*w2-25*wl3-15*»17+ll*w26+3*w34-5*w221 
+3*w442-75)*r2*r3+(174*w2+28*wl3+4*wl7+14*w26+22*w34-4*w221+14*w442-44)*r2+(-35*w2-20*wl7 
+5*w26+5*w34+5*w442-20)*r3-620*w2+160*wl3+280*wl7-60*w26-180*w34+40*w221-20*w442 
+560)*tl*t4+((73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32)*r2*r3+(20*w2 
+42*wl3+26*wl7-12*w26-16*w34+10*w221+130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 
-5*w442-55)*r3+(340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620))*tl+((6*wl3-2*wl7 
+10*w221+170)*r2*r3+(-68*wl3-44*wl7+20*w221+380)*r2+(-220*w2-60*w26+40*w34)*r3+(660*w2 
-60*w26+460*w34+60*w442) ) *t3*t4+ ( (75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442 
-7)*r2*r3+(-130*w2-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402)*r2+(-35*w2-90*wl3 
-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34 
+60*w221+10*w442+60)*t3+((128*w2-6*wl3+2*wl7-20*w26-20*w34-10*w221+8*w442-170)*r2*r3 
+(666*w2+68*wl3+44*wl7-6*w26-18*w34-20*w221+46*w442-380)*r2+(220*w2-40*wl3+80*wl7+60*w26 
-40*w34-280) *r3-660*w2+100*wl3+580*wl7+60*w26-460*w34+60*w221-60*w442+380) *t4+ (36*w2 
-26*wl3-18*wl7+10*w34-2*w221+6*w442-42)*r2*r3+(74*w2+62*wl3+26*wl7-54*w26-2*w34+26*w221 
+14*w442+286)*r2+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7 
-110*w26+170*w34-20*w221+30*w442+1060; 

cl5:=(((ll/26*wl3+l/2*wl7+3/26*w221+l/2)*r2*r3+(-10*w2-20/13*w26-10/13*w442)*r2 

-40/13*wl3*r3)*tl*t2+((43*w2+27*wl3+21*wl7+ll*w26+3*w34+3*w221+3*w442+53)*r2*r3+(174*w2 
+24*wl3-8*wl7+14*w26+22*w34+8*w221+14*w442+88)*r2+(-55*w2-10*wl3+10*wl7+25*w26+5*w34 
-10*w221+5*w442+90) *r3+ (620*w2-40*wl3+80*wl7+60*w26+180*w34+20*w442+680) ) *t 1+ ( (128*w2 
+56*wl3+28*wl7-20*w26-20*w34+4*w221+8*w442+144)*r2*r3+(666*w2+64*wl3-48*wl7-6*w26-18*w34 
-16*w221+46*w442-96)*r2+(130*w2-80*wl3+90*w26-10*w34-10*w442-160)*r3+(1060*w2-220*wl3 
+340*wl7+260*w26+140*w34+60*w221-20*w442+940)))*I+((-ll/26*wl3-l/2*wl7-3/26*w221 
-l/2)*r2*r3+(-10*w2-20/13*w26-10/13*w442)*r2)*tl*t2+((32*w2-19*wl3-17*wl7+14*w26+12*w34 
-3*w221+2*w442-73) *r2*r3+ (-130*w2+24*wl3+32*wl7-42*w26-26*w34-10*w442-40) *r2+ (-85*w2 
+20*wl7-25*w26+15*w34-5*w442+20) *r3+ (80*w2-160*wl3-280*wl7+40*w26+40*w34-40*w221-560) ) *tl 
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+(42*w2-20*wl7+26*w26+18*w34-12*w221+2*w442-72)*r2*r3+(-286*w2+108*wl3+4*wl7-62*w26 

-26*w34-28*w221-26*w442-148)*r2+(-470*w2+40*wl3-80*wl7+10*w26+30*w34-10*w442+280)*r3 

-240*w2-100*wl3-580*wl7+40*w26+400*w34-60*w221-40*w442-380; 

cl6 : = ( ( (-10*w2+4*w26+2*w34) *r2*r3+ (24*w2+24*w26+12*w34+4*w442) *r2+ (-50*wl3+10*wl7-10*w221 
-30) *r3-20*wl3-140*wl7-20*w221-540) *tl*t3*t4+( (53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221 
+3/2*w442+ll)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3 
+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 
-10*w442-10)*tl*t3+((10*w2+27*wl3+21*wl7-4*w26-2*w34+3*w221+53)*r2*r3+(-24*w2+24*wl3 
-8*wl7-24*w26-12*w34+8*w221-4*w442+88)*r2+(-85*w2+50*wl3-10*wl7-25*w26+15*w34+10*w221 
-5*w442+30)*r3+(80*w2+20*wl3+140*wl7+40*w26+40*w34+20*w221+540))*tl*t4+((-53/2*w2-ll*wl3 
-3*wl7-27/2*w26-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34 
-14*w221-4*w442-174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 
-80*w26-140*w34-20*w442+80)*tl+((36*w2+28*w26+4*w34-4*w442)*r2*r3+(-122*w2+86*w26-22*w34 
-22*w442)*r2+(-80*wl3+40*wl7-600)*r3-220*wl3+260*wl7-20*w221-1300)*t3*t4+((46*w2-25*wl3 
-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 
+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 
+10*w26-430*w34-20*w221-10*w442+460)*t3+((-36*w2+56*wl3+28*wl7-28*w26-4*w34+4*w221+4*w442 
+144)*r2*r3+(122*w2+64*wl3-48*wl7-86*w26+22*w34-16*w221+22*w442-96)*r2+(-470*w2+80*wl3 
-40*wl7+10*w26+30*w34-10*w442+600)*r3-240*w2+220*wl3-260*wl7+40*w26+400*w34+20*w221 
-40*w442+1300)*t4+((-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442-128)*r2*r3+(48*w2 
+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 
-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (3*wl3 
+9*wl7-w221-ll)*r2*r3+(-56*wl3-48*wl7-24*w221-304)*r2+(-55*w2+5*w26-15*w34+5*w442)*r3 
-60*w2+100*w26+100*w34+20*w442)*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26-31/2*w34-w221 
-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442+132)*r2+(55*w2 
-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 
-10*w221-20*w442-10)*tl*t3+((32*w2-3*wl3-9*wl7+14*w26+12*w34+w221+2*w442+ll)*r2*r3 
+(-130*w2+56*wl3+48*wl7-42*w26-26*w34+24*w221-10*w442+304)*r2+(55*w2+10*wl3-10*wl7-5*w26 
+15*w34+10*w221-5*w442-90)*r3+(60*w2+40*wl3-80*wl7-100*w26-100*w34-20*w442-680))*tl*t4 
+((73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32)*r2*r3+(20*w2+42*wl3+26*wl7 
-12*w26-16*w34+10*w221+130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442+55)*r3 
-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620)*tl+((46*wl3+38*wl7-6*w221-86)*r2*r3 
+ (-56*wl3-8*wl7-72*w221-952) *r2+ (-60*w2+20*w26+40*w34) *r3-280*w2+160*w26+120*w34) *t3*t4 
+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26 
+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3 
-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3+((42*w2-46*wl3-38*wl7 
+26*w26+18*w34+6*w221+2*w442+86)*r2*r3+(-286*w2+56*wl3+8*wl7-62*w26-26*w34+72*w221 
-26*w442+952)*r2+(60*w2+80*wl3-20*w26-40*w34+160)*r3+(280*w2+220*wl3-340*wl7-160*w26 
-120*w34-60*w221-940))*t4+(36*w2-26*wl3-18*wl7+10*w34-2*w221+6*w442-42)*r2*r3+(74*w2 
+62*wl3+26*wl7-54*w26-2*w34+26*w221+14*w442+286)*r2+(80*w2-90*wl3+10*wl7+40*w26+10*w221 
-130)*r3-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442-1060; 

phi: = [ cOl, c02, c03, c04, c05, c06, c07, c08, c09, clO, ell, cl2, cl3, cl4, cl5, cl6] ; 



166 

w2:=Sqrt(2) ; 
wl3:=Sqrt(13); 
wl7:=Sqrt(17) ; 

I:=Sqrt(-l) ; 

w26:=w2*wl3; 
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w34:=w2*wl7; 
w221:=wl3*wl7; 
w442:=w2*wl3*wl7; 

r2:=Sqrt(wl3*wl7-ll) ; 
r3:=Sqrt(15+wl7); 

tl:=Sqrt(15+(4-wl7)*r3-3*wl7) ; 

t2 : =Sqrt ( ( ( (-5*wl7+3) *wl3+ (39*wl7-65) ) *r3+ ( ( 16*wl7-72) *wl3+936) ) *t l-208*wl3+2288) ; 
t3:=Sqrt(2-w2) ; 
t4 : =Sqrt (2+t3) ; 



c01:=(-20/13*w26*r3*tl*t2+((-90*w2+60*wl3-20*wl7+30*w26+10*w34+20*w221+10*w442-60)*r3 

+(60*wl3+60*wl7+20*w221-140))*tl*t3+((350*w2+160*wl3-40*wl7+150*w26-50*w34-10*w442+760)*r3 

+(400*w2+440*wl3-600*wl7+320*w26-320*w34-40*w221-80*w442+360))*t3)*I+(-40/13*wl3 

-20/13*w26)*r3*tl*t2+((140*w2+20*w26+160)*r3-20*w2-260*wl3-500*wl7-140*w26-140*w34-60*w221 

-20*w442-1260)*tl*t3+((530*w2-60*wl3-100*wl7-30*w26-70*w34+20*w221+10*w442+620)*r3+(520*w2 

-400*wl3-1120*wl7-200*w26-520*w34+40*w442-1200))*t3; 

c02 : =( ( (10*w2-4*w26-2*w34) *r2*r3+(-24*w2-24*w26-12*w34-4*w442) *r2+ (-50*wl3+10*wl7-10*w221 

-30)*r3-20*wl3-140*wl7-20*w221-540)*tl*t3*t4+((-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221 

-3/2*w442-ll)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((-10*w2-27*wl3-21*wl7+4*w26+2*w34-3*w221-53)*r2*r3+(24*w2-24*wl3 

+8*wl7+24*w26+12*w34-8*w221+4*w442-88)*r2+(-85*w2+50*wl3-10*wl7-25*w26+15*w34+10*w221 

-5*w442+30)*r3+(80*w2+20*wl3+140*wl7+40*w26+40*w34+20*w221+540))*tl*t4+((53/2*w2+ll*wl3 

+3*wl7+27/2*w26+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34 

+14*w221+4*w442+174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((-36*w2-28*w26-4*w34+4*w442)*r2*r3+(122*w2-86*w26+22*w34 

+22*w442)*r2+(-80*wl3+40*wl7-600)*r3-220*wl3+260*wl7-20*w221-1300)*t3*t4+((-46*w2+25*wl3 

+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34-18*w221 

-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((36*w2-56*wl3-28*wl7+28*w26+4*w34-4*w221-4*w442 

-144)*r2*r3+(-122*w2-64*wl3+48*wl7+86*w26-22*w34+16*w221-22*w442+96)*r2+(-470*w2+80*wl3 

-40*wl7+10*w26+30*w34-10*w442+600)*r3-240*w2+220*wl3-260*wl7+40*w26+400*w34+20*w221 

-40*w442+1300) *t4+( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3+(-48*w2 

-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (-3*wl3 

-9*wl7+w221+ll)*r2*r3+(56*wl3+48*wl7+24*w221+304)*r2+(-55*w2+5*w26-15*w34+5*w442)*r3 

-60*w2+100*w26+100*w34+20*w442)*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7+37/2*w26+31/2*w34+w221 

+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221-3*w442-132)*r2+(55*w2 

-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((-32*w2+3*wl3+9*wl7-14*w26-12*w34-w221-2*w442-ll)*r2*r3 

+(130*w2-56*wl3-48*wl7+42*w26+26*w34-24*w221+10*w442-304)*r2+(55*w2+10*wl3-10*wl7-5*w26 

+15*w34+10*w221-5*w442-90)*r3+(60*w2+40*wl3-80*wl7-100*w26-100*w34-20*w442-680))*tl*t4 

+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2-42*wl3 

-26*wl7+12*w26+16*w34-10*w221-130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442 

+55)*r3-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620)*tl+((-46*wl3-38*wl7+6*w221 

+86) *r2*r3+ (56*wl3+8*wl7+72*w221+952) *r2+ (-60*w2+20*w26+40*w34) *r3-280*w2+160*w26 

+120*w34)*t3*t4+((75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442-7)*r2*r3+(-130*w2 

-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34 

+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3 

+((-42*w2+46*wl3+38*wl7-26*w26-18*w34-6*w221-2*w442-86)*r2*r3+(286*w2-56*wl3-8*wl7+62*w26 

+26*w34-72*w221+26*w442-952) *r2+ (60*w2+80*wl3-20*w26-40*w34+160) *r3+ (280*w2+220*wl3 
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-340*wl7-160*w26-120*w34-60*w221-940))*t4+(-36*w2+26*wl3+18*wl7-10*w34+2*w221-6*w442 
+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221-14*w442-286)*r2+(80*w2-90*wl3 
+ 10*wl7+40*w26+10*w221-130)*r3-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442 
-1060; 

c03:=(((-ll/26*wl3-l/2*wl7-3/26*w221-l/2)*r2*r3+(10*w2+20/13*w26+10/13*w442)*r2 

-40/13*wl3*r3)*tl*t2+((-43*w2-27*wl3-21*wl7-ll*w26-3*w34-3*w221-3*w442-53)*r2*r3+(-174*w2 

-24*wl3+8*wl7-14*w26-22*w34-8*w221-14*w442-88)*r2+(-55*w2-10*wl3+10*wl7+25*w26+5*w34 

-10*w221+5*w442+90) *r3+ (620*w2-40*wl3+80*wl7+60*w26+180*w34+20*w442+680) ) *t 1+ ( (-128*w2 

-56*wl3-28*wl7+20*w26+20*w34-4*w221-8*w442-144)*r2*r3+(-666*w2-64*wl3+48*wl7+6*w26+18*w34 

+16*w221-46*w442+96)*r2+(130*w2-80*wl3+90*w26-10*w34-10*w442-160)*r3+(1060*w2-220*wl3 

+340*wl7+260*w26+140*w34+60*w221-20*w442+940) ) ) *I+ ( (H/26*wl3+l/2*wl7+3/26*w221+l/2) *r2*r3 

+(10*w2+20/13*w26+10/13*w442)*r2)*tl*t2+((-32*w2+19*wl3+17*wl7-14*w26-12*w34+3*w221-2*w442 

+73)*r2*r3+(130*w2-24*wl3-32*wl7+42*w26+26*w34+10*w442+40)*r2+(-85*w2+20*wl7-25*w26 

+15*w34-5*w442+20)*r3+(80*w2-160*»13-280*wl7+40*w26+40*w34-40*w221-560))*tl+(-42*w2 

+20*wl7-26*w26-18*w34+12*w221-2*w442+72)*r2*r3+(286*w2-108*wl3-4*wl7+62*w26+26*w34 

+28*w221+26*w442+148) *r2+ (-470*w2+40*wl3-80*wl7+10*w26+30*w34-10*w442+280) *r3-240*w2 

-100*wl3-580*wl7+40*w26+400*w34-60*w221-40*w442-380; 

c04 : = ( ( (-63*w2-23*w26-19*w34-3*w442) *r2*r3+ (-64*w2+20*w34-4*w442) *r2+ (-20*wl7+140) *r3 

-100*wl3-220*wl7-20*w221-700)*tl*t3*t4+((53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221+3/2*w442 

+1 1) *r2*r3+ (75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22) *r2+ (25*w2+25*wl3+5*wl7 

+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221-10*w442 

-10)*tl*t3+((63*w2-19*wl3-17*wl7+23*w26+19*w34-3*w221+3*w442-73)*r2*r3+(64*w2+24*wl3 

+32*wl7-20*w34+4*w442-40)*r2+(55*w2+20*wl7-25*w26-5*w34-5*w442-140)*r3-620*w2+100*wl3 

+220*wl7-60*w26-180*w34+20*w221-20*w442+700)*tl*t4+((53/2*w2+ll*wl3+3*wl7+27/2*w26 

+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34+14*w221+4*w442 

+174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26+140*w34 

+20*w442-80))*tl+((-108*w2-28*w26-24*w34-8*w442)*r2*r3+(-26*w2+22*w26+26*w34-6*w442)*r2 

+(-100*wl3-20*wl7+20*w221+340)*r3-300*wl3-540*wl7+60*w221-820)*t3*t4+((46*w2-25*wl3 

-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 

+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((108*w2-20*wl7+28*w26+24*w34-12*w221+8*w442 

-72)*r2*r3+(26*w2+108*wl3+4*wl7-22*w26-26*w34-28*w221+6*w442-148)*r2+(-130*w2+100*wl3 

+20*wl7-90*w26+10*w34-20*w221+10*w442-340)*r3-1060*w2+300*wl3+540*wl7-260*w26-140*w34 

-60*w221+20*w442+820)*t4+((72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128)*r2*r3 

+(-48*w2-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(-140*w2-10*wl3-30*wl7-20*w26 

+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221+30*w442+240)))*I 

+ ( (-25*wl3-15*wl7-5*w221-75) *r2*r3+ (28*wl3+4*wl7-4*w221-44) *r2+ (35*w2-5*w26-5*w34 

-5*w442)*r3+(620*w2+60*w26+180*w34+20*w442))*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26 

-31/2*w34-w221-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442 

+132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7 

-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+((-43*w2+25*wl3+15*wl7-ll*w26-3*w34+5*w221 

-3*w442+75)*r2*r3+(-174*w2-28*wl3-4*wl7-14*w26-22*w34+4*w221-14*w442+44)*r2+(-35*w2 

-20*wl7+5*w26+5*w34+5*w442-20)*r3-620*w2+160*wl3+280*wl7-60*w26-180*w34+40*w221-20*w442 

+560)*tl*t4+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2 

-42*wl3-26*wl7+12*w26+16*w34-10*w221-130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 

-5*w442-55)*r3+(340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620))*tl+((-6*wl3+2*wl7 

-10*w221-170)*r2*r3+(68*wl3+44*wl7-20*w221-380)*r2+(-220*w2-60*w26+40*w34)*r3+(660*w2 

-60*w26+460*w34+60*w442))*t3*t4+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442 

+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3 

-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34 

+60*w221+10*w442+60)*t3+((-128*w2+6*wl3-2*wl7+20*w26+20*w34+10*w221-8*w442+170)*r2*r3 

+ (-666*w2-68*wl3-44*wl7+6*w26+18*w34+20*w221-46*w442+380) *r2+ (220*w2-40*wl3+80*wl7+60*w26 

-40*w34-280) *r3-660*w2+100*wl3+580*wl7+60*w26-460*w34+60*w221-60*w442+380) *t4+ (-36*w2 
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+26*wl3+18*wl7-10*w34+2*w221-6*w442+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221 

-14*w442-286)*r2+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7 

-110*w26+170*w34-20*w221+30*w442+1060; 

c05:=(-20/13*w26*r3*tl*t2+((20*w2+30*wl3+30*wl7+12*w26+16*w34+2*w221+42)*r2*r3+(-110*w2 

-32*wl3-56*wl7-14*w26-42*w34-16*w221-10*w442-216)*r2)*tl*t3+((-20*w2+102*wl3+66*wl7+48*w26 

+44*w34-2*w221+58)*r2*r3+(-640*w2+8*wl3-56*wl7-16*w26-8*w34-88*w221-40*w442 

-1048) *r2) *t3) *I+(40/13*wl3-20/13*w26) *r3*tl*t2+( (42*w2+44*wl3+32*wl7+10*w26+10*w34+8*w221 

+2*w442+148)*r2*r3+(-154*w2-52*wl3-36*wl7-66*w26-38*w34+4*w221-14*w442+84)*r2)*tl*t3 

+((6*w2+6*wl3+18*wl7-2*w26+14*w34+22*w221+6*w442+242)*r2*r3+(-164*w2-176*wl3-48*wl7 

-148*w26-4*w34+48*w221-4*w442+528) *r2) *t3 ; 

c06 : = ( ( (-63*w2-23*w26-19*w34-3*w442) *r2*r3+ (-64*w2+20*w34-4*w442) *r2+ (20*wl7-140) *r3 

+(100*wl3+220*wl7+20*w221+700))*tl*t3*t4+((-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221 

-3/2*w442-ll)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*»17-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((63*w2-19*wl3-17*wl7+23*w26+19*w34-3*w221+3*w442-73)*r2*r3+(64*w2 

+24*wl3+32*wl7-20*w34+4*w442-40) *r2+ (-55*w2-20*wl7+25*w26+5*w34+5*w442+140) *r3+ (620*w2 

-100*wl3-220*wl7+60*w26+180*w34-20*w221+20*w442-700))*tl*t4+((-53/2*w2-ll*wl3-3*wl7 

-27/2*w26-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34-14*w221 

-4*w442-174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26 

+140*w34+20*w442-80) ) *tl+ ( (-108*w2-28*w26-24*w34-8*w442) *r2*r3+ (-26*w2+22*w26+26*w34 

-6*w442) *r2+ (100*wl3+20*wl7-20*w221-340) *r3+ (300*wl3+540*wl7-60*w221+820) ) *t3*t4+ ( (-46*w2 

+25*wl3+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34 

-18*w221-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3 

-420*wl7+10*w26-430*w34-20*w221-10*w442+460)*t3+((108*w2-20*wl7+28*w26+24*w34-12*w221 

+8*w442-72)*r2*r3+(26*w2+108*wl3+4*wl7-22*w26-26*w34-28*w221+6*w442-148)*r2+(130*w2 

-100*wl3-20*wl7+90*w26-10*w34+20*w221-10*w442+340)*r3+(1060*w2-300*wl3-540*wl7+260*w26 

+140*w34+60*w221-20*w442-820))*t4+((-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442 

-128)*r2*r3+(48*w2+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(-140*w2-10*wl3 

-30*wl7-20*w26+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221 

+30*w442+240) ) ) *I+ ( (-25*wl3-15*wl7-5*w221-75) *r2*r3+ (28*wl3+4*wl7-4*w221-44) *r2+ (-35*w2 

+5*w26+5*w34+5*w442)*r3-620*w2-60*w26-180*w34-20*w442)*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7 

+37/2*w26+31/2*w34+w221+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221 

-3*w442-132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3 

-70*wl7-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+((-43*w2+25*wl3+15*wl7-ll*w26-3*w34 

+5*w221-3*w442+75)*r2*r3+(-174*w2-28*wl3-4*wl7-14*w26-22*w34+4*w221-14*w442+44)*r2+(35*w2 

+20*wl7-5*w26-5*w34-5*w442+20)*r3+(620*w2-160*wl3-280*wl7+60*w26+180*w34-40*w221+20*w442 

-560))*tl*t4+((73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32)*r2*r3+(20*w2 

+42*wl3+26*wl7-12*w26-16*w34+10*w221+130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 

-5*w442-55)*r3+(340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620))*tl+((-6*wl3+2*wl7 

-10*w221-170)*r2*r3+(68*wl3+44*wl7-20*w221-380)*r2+(220*w2+60*w26-40*w34)*r3-660*w2+60*w26 

-460*w34-60*w442) *t3*t4+ ( (75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442-7) *r2*r3 

+ (-130*w2-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402) *r2+ (-35*w2-90*wl3-10*wl7-55*w26 

-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442 

+60)*t3+((-128*w2+6*wl3-2*wl7+20*w26+20*w34+10*w221-8*w442+170)*r2*r3+(-666*w2-68*wl3 

-44*wl7+6*w26+18*w34+20*w221-46*w442+380) *r2+ (-220*w2+40*wl3-80*wl7-60*w26+40*w34+280) *r3 

+(660*w2-100*wl3-580*wl7-60*w26+460*w34-60*w221+60*w442-380))*t4+(36*w2-26*wl3-18*wl7 

+10*w34-2*w221+6*w442-42) *r2*r3+ (74*w2+62*wl3+26*wl7-54*w26-2*w34+26*w221+14*w442+286) *r2 

+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7-110*w26+170*w34 

-20*w221+30*w442+1060 ; 

c07:=(((ll/26*wl3+l/2*wl7+3/26*w221+l/2)*r2*r3+(-10*w2-20/13*w26-10/13*w442)*r2 

-40/13*wl3*r3)*tl*t2+((-43*w2-27*wl3-21*wl7-ll*w26-3*w34-3*w221-3*w442-53)*r2*r3+(-174*w2 
-24*wl3+8*wl7-14*w26-22*w34-8*w221-14*w442-88)*r2+(55*w2+10*wl3-10*wl7-25*w26-5*w34 
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+10*w221-5*w442-90) *r3-620*w2+40*wl3-80*wl7-60*w26-180*w34-20*w442-680) *tl+ ( (-128*w2 

-56*wl3-28*wl7+20*w26+20*w34-4*w221-8*w442-144)*r2*r3+(-666*w2-64*wl3+48*wl7+6*w26+18*w34 

+16*w221-46*w442+96)*r2+(-130*w2+80*wl3-90*w26+10*w34+10*w442+160)*r3-1060*w2+220*wl3 

-340*wl7-260*w26-140*w34-60*w221+20*w442-940))*I+((-ll/26*wl3-l/2*wl7-3/26*w221-l/2)*r2*r3 

+(-10*w2-20/13*w26-10/13*w442)*r2)*tl*t2+((-32*w2+19*wl3+17*wl7-14*w26-12*w34+3*w221-2*w442 

+73)*r2*r3+(130*w2-24*wl3-32*wl7+42*w26+26*w34+10*w442+40)*r2+(85*w2-20*wl7+25*w26-15*w34 

+5*w442-20)*r3-80*w2+160*wl3+280*wl7-40*w26-40*w34+40*w221+560)*tl+(-42*w2+20*wl7-26*w26 

-18*w34+12*w221-2*w442+72)*r2*r3+(286*w2-108*wl3-4*wl7+62*w26+26*w34+28*w221+26*w442 

+148) *r2+ (470*w2-40*wl3+80*wl7-10*w26-30*w34+10*w442-280) *r3+240*w2+100*wl3+580*wl7 

-40*w26-400*w34+60*w22 l+40*w442+380 ; 

c08:=(((10*w2-4*w26-2*w34)*r2*r3+(-24*w2-24*w26-12*w34-4*w442)*r2+(50*wl3-10*wl7+10*w221 

+30)*r3+(20*wl3+140*wl7+20*w221+540))*tl*t3*t4+((53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221 

+3/2*w442+ll)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((-10*w2-27*wl3-21*wl7+4*w26+2*w34-3*w221-53)*r2*r3+(24*w2-24*wl3 

+8*wl7+24*w26+12*w34-8*w221+4*w442-88)*r2+(85*w2-50*wl3+10*wl7+25*w26-15*w34-10*w221 

+5*w442-30)*r3-80*w2-20*wl3-140*wl7-40*w26-40*w34-20*w221-540)*tl*t4+((-53/2*w2-ll*wl3 

-3*wl7-27/2*w26-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34 

-14*w221-4*w442-174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((-36*w2-28*w26-4*w34+4*w442)*r2*r3+(122*w2-86*w26+22*w34 

+22*w442) *r2+ (80*wl3-40*wl7+600) *r3+ (220*wl3-260*wl7+20*w221+1300) ) *t3*t4+ ( (46*w2-25*wl3 

-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 

+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((36*w2-56*wl3-28*wl7+28*w26+4*w34-4*w221-4*w442 

-144) *r2*r3+ (-122*w2-64*wl3+48*wl7+86*w26-22*w34+16*w221-22*w442+96) *r2+ (470*w2-80*wl3 

+40*wl7-10*w26-30*w34+10*w442-600)*r3+(240*w2-220*wl3+260*wl7-40*w26-400*w34-20*w221 

+40*w442-1300))*t4+((-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442-128)*r2*r3+(48*w2 

+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (-3*wl3 

-9*wl7+w221+ll)*r2*r3+(56*wl3+48*wl7+24*w221+304)*r2+(55*w2-5*w26+15*w34-5*w442)*r3 

+(60*w2-100*w26-100*w34-20*w442))*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26-31/2*w34-w221 

-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442+132)*r2+(55*w2 

-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((-32*w2+3*wl3+9*wl7-14*w26-12*w34-w221-2*w442-ll)*r2*r3 

+(130*w2-56*wl3-48*wl7+42*w26+26*w34-24*w221+10*w442-304)*r2+(-55*w2-10*wl3+10*wl7+5*w26 

-15*w34-10*w221+5*w442+90)*r3-60*w2-40*wl3+80*wl7+100*w26+100*w34+20*w442+680)*tl*t4 

+ ( (73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32) *r2*r3+ (20*w2+42*wl3+26*wl7 

-12*w26-16*w34+10*w221+130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442+55)*r3 

-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620) *tl+( (-46*wl3-38*wl7+6*w221+86) *r2*r3 

+(56*wl3+8*wl7+72*w221+952) *r2+(60*w2-20*w26-40*w34) *r3+(280*w2-160*w26-120*w34) ) *t3*t4 

+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26 

+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3 

-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3+((-42*w2+46*wl3+38*wl7 

-26*w26-18*w34-6*w221-2*w442-86)*r2*r3+(286*w2-56*wl3-8*wl7+62*w26+26*w34-72*w221+26*w442 

-952)*r2+(-60*w2-80*wl3+20*w26+40*w34-160)*r3-280*w2-220*wl3+340*wl7+160*w26+120*w34 

+60*w221+940)*t4+(36*w2-26*wl3-18*wl7+10*w34-2*w221+6*w442-42)*r2*r3+(74*w2+62*wl3 

+26*wl7-54*w26-2*w34+26*w221+14*w442+286)*r2+(80*w2-90*wl3+10*wl7+40*w26+10*w221-130)*r3 

-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442-1060; 

c09:=(-20/13*w26*r3*tl*t2+((90*w2-60*wl3+20*wl7-30*w26-10*w34-20*w221-10*w442+60)*r3-60*wl3 

-60*wl7-20*w221+140)*tl*t3+((-350*w2-160*wl3+40*wl7-150*w26+50*w34+10*w442-760)*r3-400*w2 
-440*wl3+600*wl7-320*w26+320*w34+40*w221+80*w442-360)*t3)*I+(-40/13*wl3-20/13*w26)*r3*tl*t2 
+((-140*w2-20*w26-160)*r3+(20*w2+260*wl3+500*wl7+140*w26+140*w34+60*w221+20*w442 
+1260))*tl*t3+((-530*w2+60*wl3+100*wl7+30*w26+70*w34-20*w221-10*w442-620)*r3-520*w2 
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+400*wl3+1120*wl7+200*w26+520*w34-40*w442+1200) *t3 ; 

clO : = ( ( (-10*w2+4*w26+2*w34) *r2*r3+ (24*w2+24*w26+12*w34+4*w442) *r2+ (50*wl3-10*wl7+10*w221 

+30) *r3+ (20*wl3+140*wl7+20*w221+540) ) *t I*t3*t4+ ( (-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221 

-3/2*w442-ll)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3 

+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 

-10*w442-10)*tl*t3+((10*w2+27*wl3+21*wl7-4*w26-2*w34+3*w221+53)*r2*r3+(-24*w2+24*wl3 

-8*wl7-24*w26-12*w34+8*w221-4*w442+88)*r2+(85*w2-50*wl3+10*wl7+25*w26-15*w34-10*w221 

+5*w442-30) *r3-80*w2-20*wl3-140*wl7-40*w26-40*w34-20*w221-540) *t l*t4+ ( (53/2*w2+ll*wl3 

+3*wl7+27/2*w26+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34 

+14*w221+4*w442+174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((36*w2+28*w26+4*w34-4*w442)*r2*r3+(-122*w2+86*w26-22*w34 

-22*w442) *r2+ (80*wl3-40*wl7+600) *r3+ (220*wl3-260*wl7+20*w221+1300) ) *t3*t4+ ( (-46*w2+25*wl3 

+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34-18*w221 

-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((-36*w2+56*wl3+28*wl7-28*w26-4*w34+4*w221+4*w442 

+144)*r2*r3+(122*w2+64*wl3-48*wl7-86*w26+22*w34-16*w221+22*w442-96)*r2+(470*w2-80*wl3 

+40*wl7-10*w26-30*w34+10*w442-600)*r3+(240*w2-220*wl3+260*wl7-40*w26-400*w34-20*w221 

+40*w442-1300) ) *t4+ ( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3+ (-48*w2 

-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470) *r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240) ) *I+ ( (3*wl3 

+9*wl7-w221-ll)*r2*r3+(-56*wl3-48*wl7-24*w221-304)*r2+(55*w2-5*»26+15*w34-5*w442)*r3 

+(60*w2-100*w26-100*w34-20*w442))*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7+37/2*w26+31/2*w34+w221 

+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221-3*w442-132)*r2+(55*w2 

-5*wl3-5*»17-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((32*w2-3*wl3-9*wl7+14*w26+12*w34+w221+2*w442+ll)*r2*r3 

+(-130*w2+56*wl3+48*wl7-42*w26-26*w34+24*w221-10*w442+304)*r2+(-55*w2-10*wl3+10*wl7+5*w26 

-15*w34-10*w221+5*w442+90)*r3-60*w2-40*wl3+80*wl7+100*w26+100*w34+20*w442+680)*tl*t4 

+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2-42*wl3 

-26*wl7+12*w26+16*w34-10*w221-130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442 

+55)*r3-340*w2-60*wl3-180*wl7+20*w26-40*w34-20*w221-620)*tl+((46*wl3+38*wl7-6*w221 

-86)*r2*r3+(-56*wl3-8*wl7-72*w221-952)*r2+(60*w2-20*w26-40*w34)*r3+(280*w2-160*w26 

-120*w34))*t3*t4+((75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442-7)*r2*r3+(-130*w2 

-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34 

+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3 

+ ( (42*w2-46*wl3-38*wl7+26*w26+18*w34+6*w221+2*w442+86) *r2*r3+ (-286*w2+56*wl3+8*wl7-62*w26 

-26*w34+72*w221-26*w442+952)*r2+(-60*w2-80*wl3+20*w26+40*w34-160)*r3-280*w2-220*wl3 

+340*wl7+160*w26+120*w34+60*w221+940)*t4+(-36*w2+26*wl3+18*wl7-10*w34+2*w221-6*w442 

+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221-14*w442-286)*r2+(80*w2-90*wl3 

+ 10*wl7+40*w26+10*w221-130)*r3-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442 

-1060; 

cll:=(((-ll/26*wl3-l/2*wl7-3/26*w221-l/2)*r2*r3+(10*w2+20/13*w26+10/13*w442)*r2 

-40/13*wl3*r3)*tl*t2+((43*w2+27*wl3+21*wl7+ll*w26+3*w34+3*w221+3*w442+53)*r2*r3+(174*w2 

+24*wl3-8*wl7+14*w26+22*w34+8*w221+14*w442+88)*r2+(55*w2+10*wl3-10*wl7-25*w26-5*w34 

+10*w221-5*w442-90) *r3-620*w2+40*wl3-80*wl7-60*w26-180*w34-20*w442-680) *t 1+ ( ( 128*w2 

+56*wl3+28*wl7-20*w26-20*w34+4*w221+8*w442+144)*r2*r3+(666*w2+64*wl3-48*wl7-6*w26-18*w34 

-16*w221+46*w442-96)*r2+(-130*w2+80*wl3-90*w26+10*w34+10*w442+160)*r3-1060*w2+220*wl3 

-340*wl7-260*w26-140*w34-60*w221+20*w442-940))*I+((ll/26*wl3+l/2*wl7+3/26*w221+l/2)*r2*r3 

+(10*w2+20/13*w26+10/13*w442)*r2)*tl*t2+((32*w2-19*wl3-17*wl7+14*w26+12*w34-3*w221+2*w442 

-73) *r2*r3+ (-130*w2+24*wl3+32*wl7-42*w26-26*w34-10*w442-40) *r2+ (85*w2-20*wl7+25*w26 

-15*w34+5*w442-20)*r3-80*w2+160*wl3+280*wl7-40*w26-40*w34+40*w221+560)*tl+(42*w2-20*wl7 

+26*w26+18*w34-12*w221+2*w442-72)*r2*r3+(-286*w2+108*wl3+4*wl7-62*w26-26*w34-28*w221 

-26*w442-148) *r2+ (470*w2-40*wl3+80*wl7-10*w26-30*w34+10*w442-280) *r3+240*w2+100*wl3 

+580*wl7-40*w26-400*w34+60*w221+40*w442+380; 
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cl2 : = ( ( (63*w2+23*w26+19*w34+3*w442) *r2*r3+ (64*w2-20*w34+4*w442) *r2+ (20*wl7-140) *r3+ (100*wl3 
+220*wl7+20*w221+700))*tl*t3*t4+((53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221+3/2*w442 
+Il)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3+5*wl7 
+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221-10*w442 
-10)*tl*t3+((-63*w2+19*wl3+17*wl7-23*w26-19*w34+3*w221-3*w442+73)*r2*r3+(-64*w2-24*wl3 
-32*wl7+20*w34-4*w442+40)*r2+(-55*w2-20*wl7+25*w26+5*w34+5*w442+140)*r3+(620*w2-100*wl3 
-220*wl7+60*w26+180*w34-20*w221+20*w442-700))*tl*t4+((53/2*w2+ll*wl3+3*wl7+27/2*w26 
+21/2*w34+3*w221+3/2*w442+43)*r2*r3+(44*w2+14*wl3+22*wl7+12*w26-4*w34+14*w221+4*w442 
+174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26+140*w34 
+20*w442-80) ) *t 1+ ( ( 108*w2+28*w26+24*w34+8*w442) *r2*r3+ (26*w2-22*w26-26*w34+6*w442) *r2 
+(100*wl3+20*wl7-20*w221-340)*r3+(300*wl3+540*wl7-60*w221+820))*t3*t4+((46*w2-25*wl3 
-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 
+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 
+10*w26-430*w34-20*w221-10*w442+460)*t3+((-108*w2+20*wl7-28*w26-24*w34+12*w221-8*w442 
+72)*r2*r3+(-26*w2-108*wl3-4*wl7+22*w26+26*w34+28*w221-6*w442+148)*r2+(130*w2-100*wl3 
-20*wl7+90*w26-10*w34+20*w221-10*w442+340)*r3+(1060*w2-300*wl3-540*wl7+260*w26+140*w34 
+60*w221-20*w442-820) ) *t4+ ( (72*w2-20*wl3-20*wl7+28*w26+14*w34+8*w221+2*w442+128) *r2*r3 
+(-48*w2-6*wl3-18*wl7+32*w26-24*w34+46*w221-8*w442+666)*r2+(-140*w2-10*wl3-30*wl7-20*w26 
+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221+30*w442+240)))*I 
+ ( (25*wl3+15*wl7+5*w221+75) *r2*r3+ (-28*wl3-4*wl7+4*w221+44) *r2+ (-35*w2+5*w26+5*w34 
+5*w442)*r3-620*w2-60*w26-180*w34-20*w442)*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26 
-31/2*w34-w221-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442 
+132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7 
-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+( (43*w2-25*wl3-15*wl7+ll*w26+3*w34-5*w221 
+3*w442-75)*r2*r3+(174*w2+28*wl3+4*wl7+14*w26+22*w34-4*w221+14*w442-44)*r2+(35*w2+20*wl7 
-5*w26-5*w34-5*w442+20)*r3+(620*w2-160*wl3-280*wl7+60*w26+180*w34-40*w221+20*w442 
-560))*tl*t4+((-73/2*w2+14*wl3+12*wl7-19/2*w26-17/2*w34+2*w221-3/2*w442+32)*r2*r3+(-20*w2 
-42*wl3-26*wl7+12*w26+16*w34-10*w221-130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 
-5*w442-55) *r3+ (340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620) ) *t 1+ ( (6*wl3-2*wl7 
+10*w221+170)*r2*r3+(-68*wl3-44*wl7+20*w221+380)*r2+(220*w2+60*w26-40*w34)*r3-660*w2 
+60*w26-460*w34-60*w442) *t3*t4+ ( (-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442 
+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3 
-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34 
+60*w221+10*w442+60)*t3+((128*w2-6*wl3+2*wl7-20*w26-20*w34-10*w221+8*w442-170)*r2*r3 
+(666*w2+68*wl3+44*wl7-6*w26-18*w34-20*w221+46*w442-380)*r2+(-220*w2+40*wl3-80*wl7-60*w26 
+40*w34+280)*r3+(660*w2-100*wl3-580*wl7-60*w26+460*w34-60*w221+60*w442-380))*t4+(-36*w2 
+26*wl3+18*wl7-10*w34+2*w221-6*w442+42)*r2*r3+(-74*w2-62*wl3-26*wl7+54*w26+2*w34-26*w221 
-14*w442-286)*r2+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7 
-110*w26+170*w34-20*w221+30*w442+1060; 

cl3:=(-20/13*w26*r3*tl*t2+((-20*w2-30*wl3-30*wl7-12*w26-16*w34-2*w221-42)*r2*r3+(110*w2 

+32*wl3+56*wl7+14*w26+42*w34+16*w221+10*w442+216)*r2)*tl*t3+((20*w2-102*wl3-66*wl7-48*w26 
-44*w34+2*w221-58)*r2*r3+(640*w2-8*wl3+56*wl7+16*w26+8*w34+88*w221+40*w442 

+1048) *r2) *t3) *I+(40/13*wl3-20/13*w26) *r3*tl*t2+( (-42*w2-44*wl3-32*wl7-10*w26-10*w34-8*w221 
-2*w442-148)*r2*r3+(154*w2+52*wl3+36*wl7+66*w26+38*w34-4*w221+14*w442-84)*r2)*tl*t3 
+((-6*w2-6*wl3-18*wl7+2*w26-14*w34-22*w221-6*w442-242)*r2*r3+(164*w2+176*wl3+48*wl7 
+ 148*w26+4*w34-48*w221+4*w442-528) *r2) *t3 ; 

cl4 : = ( ( (63*w2+23*w26+19*w34+3*w442) *r2*r3+(64*w2-20*w34+4*w442) *r2+ (-20*wl7+140) *r3-100*wl3 
-220*wl7-20*w221-700)*tl*t3*t4+((-53/2*w2+wl3+3*wl7-7/2*w26-l/2*w34-w221-3/2*w442 
-H)*r2*r3+(-75*w2+26*wl3-2*wl7+5*w26-5*w34+2*w221-5*w442+22)*r2+(25*w2+25*wl3+5*wl7 
+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221-10*w442 
-10)*tl*t3+((-63*w2+19*wl3+17*wl7-23*w26-19*w34+3*w221-3*w442+73)*r2*r3+(-64*w2-24*wl3 
-32*wl7+20*w34-4*w442+40) *r2+ (55*w2+20*wl7-25*w26-5*w34-5*w442-140) *r3-620*w2+100*wl3 
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+220*wl7-60*w26-180*w34+20*w221-20*w442+700)*tl*t4+((-53/2*w2-ll*wl3-3*wl7-27/2*w26 

-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34-14*w221-4*w442 

-174)*r2+(-10*w2+25*wl3-15*wl7-10*w34+5*w221+85)*r3+(280*w2-40*wl3-40*wl7+80*w26+140*w34 

+20*w442-80) ) *t 1+ ( ( 108*w2+28*w26+24*w34+8*w442) *r2*r3+ (26*w2-22*w26-26*w34+6*w442) *r2 

+(-100*wl3-20*wl7+20*w221+340)*r3-300*wl3-540*wl7+60*w221-820)*t3*t4+((-46*w2+25*wl3 

+15*wl7+24*w26+12*w34-3*w221-6*w442-13)*r2*r3+(-394*w2+66*wl3-2*wl7+46*w26-2*w34-18*w221 

-34*w442-238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((-108*w2+20*wl7-28*w26-24*w34+12*w221-8*w442 

+72)*r2*r3+(-26*w2-108*wl3-4*wl7+22*w26+26*w34+28*w221-6*w442+148)*r2+(-130*w2+100*wl3 

+20*wl7-90*w26+10*w34-20*w221+10*w442-340)*r3-1060*w2+300*wl3+540*wl7-260*w26-140*w34 

-60*w221+20*w442+820) *t4+ ( (-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442-128) *r2*r3 

+(48*w2+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(-140*w2-10*wl3-30*wl7-20*w26 

+40*w34+10*w221+470)*r3+(190*w2-40*wl3-400*wl7+50*w26+290*w34+40*w221+30*w442+240)))*I 

+((25*wl3+15*wl7+5*w221+75)*r2*r3+(-28*wl3-4*wl7+4*w221+44)*r2+(35*w2-5*w26-5*w34 

-5*w442)*r3+(620*w2+60*w26+180*w34+20*w442))*tl*t3*t4+((95/2*w2+ll*wl3+13*wl7+37/2*w26 

+31/2*w34+w221+5/2*w442+31)*r2*r3+(-33*w2-40*wl3-40*wl7-21*w26-23*w34-12*w221-3*w442 

-132)*r2+(55*w2-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7 

-80*w26-140*w34-10*w221-20*w442-10)*tl*t3+( (43*w2-25*wl3-15*wl7+ll*w26+3*w34-5*w221 

+3*w442-75)*r2*r3+(174*w2+28*wl3+4*wl7+14*w26+22*w34-4*w221+14*w442-44)*r2+(-35*w2-20*wl7 

+5*w26+5*w34+5*w442-20)*r3-620*w2+160*wl3+280*wl7-60*w26-180*w34+40*w221-20*w442 

+560)*tl*t4+((73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32)*r2*r3+(20*w2 

+42*wl3+26*wl7-12*w26-16*w34+10*w221 + 130)*r2+(45*w2+25*wl3+5*wl7-5*w26+5*w34+5*w221 

-5*w442-55) *r3+ (340*w2+60*wl3+180*wl7-20*w26+40*w34+20*w221+620) ) *t 1+ ( (6*wl3-2*wl7 

+10*w221+170)*r2*r3+(-68*wl3-44*wl7+20*w221+380)*r2+(-220*w2-60*w26+40*w34)*r3+(660*w2 

-60*w26+460*w34+60*w442) ) *t3*t4+ ( (75*w2+23*wl3+29*wl7+27*w26+21*w34+3*w221+5*w442 

-7)*r2*r3+(-130*w2-82*wl3-6*wl7-42*w26-26*w34-22*w221-10*w442-402)*r2+(-35*w2-90*wl3 

-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3-390*w2-260*wl3-100*wl7-210*w26-130*w34 

+60*w221+10*w442+60)*t3+((128*w2-6*wl3+2*wl7-20*w26-20*w34-10*w221+8*w442-170)*r2*r3 

+(666*w2+68*wl3+44*wl7-6*w26-18*w34-20*w221+46*w442-380)*r2+(220*w2-40*wl3+80*wl7+60*w26 

-40*w34-280) *r3-660*w2+100*wl3+580*wl7+60*w26-460*w34+60*w221-60*w442+380) *t4+ (36*w2 

-26*wl3-18*wl7+10*w34-2*w221+6*w442-42)*r2*r3+(74*w2+62*wl3+26*wl7-54*w26-2*w34+26*w221 

+14*w442+286)*r2+(-80*w2+90*wl3-10*wl7-40*w26-10*w221+130)*r3+470*w2+260*wl3+140*wl7 

-110*w26+170*w34-20*w221+30*w442+1060; 

cl5 : =( ( (H/26*wl3+l/2*wl7+3/26*w221+l/2) *r2*r3+(-10*w2-20/13*w26-10/13*w442) *r2 

-40/13*wl3*r3)*tl*t2+((43*w2+27*wl3+21*wl7+ll*w26+3*w34+3*w221+3*w442+53)*r2*r3+(174*w2 

+24*wl3-8*wl7+14*w26+22*w34+8*w221+14*w442+88)*r2+(-55*w2-10*wl3+10*wl7+25*w26+5*w34 

-10*w221+5*w442+90) *r3+ (620*w2-40*wl3+80*wl7+60*w26+180*w34+20*w442+680) ) *t 1+ ( (128*w2 

+56*wl3+28*wl7-20*w26-20*w34+4*w221+8*w442+144)*r2*r3+(666*w2+64*wl3-48*wl7-6*w26-18*w34 

-16*w221+46*w442-96)*r2+(130*w2-80*wl3+90*w26-10*w34-10*w442-160)*r3+(1060*w2-220*wl3 

+340*wl7+260*w26+140*w34+60*w221-20*w442+940) ) ) *I+ ( (-11/26*w13-1/2*w17-3/26*w221 

-l/2)*r2*r3+(-10*w2-20/13*w26-10/13*w442)*r2)*tl*t2+((32*w2-19*wl3-17*wl7+14*w26+12*w34 

-3*w221+2*w442-73)*r2*r3+(-130*w2+24*wl3+32*wl7-42*w26-26*w34-10*w442-40)*r2+(-85*w2 

+20*wl7-25*w26+15*w34-5*w442+20)*r3+(80*w2-160*wl3-280*wl7+40*w26+40*w34-40*w221-560))*tl 

+(42*w2-20*wl7+26*w26+18*w34-12*w221+2*w442-72)*r2*r3+(-286*w2+108*wl3+4*wl7-62*w26 

-26*w34-28*w221-26*w442-148)*r2+(-470*w2+40*wl3-80*wl7+10*w26+30*w34-10*w442+280)*r3 

-240*w2-100*wl3-580*wl7+40*w26+400*w34-60*w221-40*w442-380; 

cl6 : = ( ( (-10*w2+4*w26+2*w34) *r2*r3+ (24*w2+24*w26+12*w34+4*w442) *r2+ (-50*wl3+10*wl7-10*w221 
-30) *r3-20*wl3-140*wl7-20*w221-540) *tl*t3*t4+( (53/2*w2-wl3-3*wl7+7/2*w26+l/2*w34+w221 
+3/2*w442+ll)*r2*r3+(75*w2-26*wl3+2*wl7-5*w26+5*w34-2*w221+5*w442-22)*r2+(25*w2+25*wl3 
+5*wl7+15*w26-5*w34+5*w221+5*w442+25)*r3-350*w2-70*wl3-70*wl7-50*w26-110*w34-10*w221 
-10*w442-10)*tl*t3+((10*w2+27*wl3+21*wl7-4*w26-2*w34+3*w221+53)*r2*r3+(-24*w2+24*wl3 
-8*wl7-24*w26-12*w34+8*w221-4*w442+88)*r2+(-85*w2+50*wl3-10*wl7-25*w26+15*w34+10*w221 
-5*w442+30)*r3+(80*w2+20*wl3+140*wl7+40*w26+40*w34+20*w221+540))*tl*t4+((-53/2*w2-ll*wl3 
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-3*wl7-27/2*w26-21/2*w34-3*w221-3/2*w442-43)*r2*r3+(-44*w2-14*wl3-22*wl7-12*w26+4*w34 

-14*w221-4*w442-174)*r2+(10*w2-25*wl3+15*wl7+10*w34-5*w221-85)*r3-280*w2+40*wl3+40*wl7 

-80*w26-140*w34-20*w442+80)*tl+((36*w2+28*w26+4*w34-4*w442)*r2*r3+(-122*w2+86*w26-22*w34 

-22*w442)*r2+(-80*wl3+40*wl7-600)*r3-220*wl3+260*wl7-20*w221-1300)*t3*t4+((46*w2-25*wl3 

-15*wl7-24*w26-12*w34+3*w221+6*w442+13)*r2*r3+(394*w2-66*wl3+2*wl7-46*w26+2*w34+18*w221 

+34*w442+238) *r2+ (345*w2+60*wl3-60*wl7+25*w26-35*w34+5*w442+440) *r3-210*w2+60*wl3-420*wl7 

+10*w26-430*w34-20*w221-10*w442+460)*t3+((-36*w2+56*wl3+28*wl7-28*w26-4*w34+4*w221+4*w442 

+144)*r2*r3+(122*w2+64*wl3-48*wl7-86*w26+22*w34-16*w221+22*w442-96)*r2+(-470*w2+80*wl3 

-40*wl7+10*w26+30*w34-10*w442+600)*r3-240*w2+220*wl3-260*wl7+40*w26+400*w34+20*w221 

-40*w442+1300)*t4+((-72*w2+20*wl3+20*wl7-28*w26-14*w34-8*w221-2*w442-128)*r2*r3+(48*w2 

+6*wl3+18*wl7-32*w26+24*w34-46*w221+8*w442-666)*r2+(140*w2+10*wl3+30*wl7+20*w26-40*w34 

-10*w221-470)*r3-190*w2+40*wl3+400*wl7-50*w26-290*w34-40*w221-30*w442-240))*I+((3*wl3 

+9*wl7-w221-ll)*r2*r3+(-56*wl3-48*wl7-24*w221-304)*r2+(-55*w2+5*w26-15*w34+5*w442)*r3 

-60*w2+100*w26+100*w34+20*w442)*tl*t3*t4+((-95/2*w2-ll*wl3-13*wl7-37/2*w26-31/2*w34-w221 

-5/2*w442-31)*r2*r3+(33*w2+40*wl3+40*wl7+21*w26+23*w34+12*w221+3*w442+132)*r2+(55*w2 

-5*wl3-5*wl7-15*w26+5*w34-5*w221-5*w442+115)*r3-280*w2-70*wl3-70*wl7-80*w26-140*w34 

-10*w221-20*w442-10)*tl*t3+((32*w2-3*wl3-9*wl7+14*w26+12*w34+w221+2*w442+ll)*r2*r3 

+(-130*w2+56*wl3+48*wl7-42*w26-26*w34+24*w221-10*w442+304)*r2+(55*w2+10*wl3-10*wl7-5*w26 

+15*w34+10*w221-5*w442-90)*r3+(60*w2+40*wl3-80*wl7-100*w26-100*w34-20*w442-680))*tl*t4 

+ ( (73/2*w2-14*wl3-12*wl7+19/2*w26+17/2*w34-2*w221+3/2*w442-32) *r2*r3+ (20*w2+42*wl3+26*wl7 

-12*w26-16*w34+10*w221+130)*r2+(-45*w2-25*wl3-5*wl7+5*w26-5*w34-5*w221+5*w442+55)*r3 

-340*w2-60*»13-180*wl7+20*w26-40*w34-20*w221-620)*tl+((46*wl3+38*wl7-6*w221-86)*r2*r3 

+ (-56*wl3-8*wl7-72*w221-952) *r2+ (-60*w2+20*w26+40*w34) *r3-280*w2+160*w26+120*w34) *t3*t4 

+((-75*w2-23*wl3-29*wl7-27*w26-21*w34-3*w221-5*w442+7)*r2*r3+(130*w2+82*wl3+6*wl7+42*w26 

+26*w34+22*w221+10*w442+402)*r2+(-35*w2-90*wl3-10*wl7-55*w26-15*w34+10*w221+5*w442+90)*r3 

-390*w2-260*wl3-100*wl7-210*w26-130*w34+60*w221+10*w442+60)*t3+((42*w2-46*wl3-38*wl7 

+26*w26+18*w34+6*w221+2*w442+86)*r2*r3+(-286*w2+56*wl3+8*wl7-62*w26-26*w34+72*w221 

-26*w442+952)*r2+(60*w2+80*wl3-20*w26-40*w34+160)*r3+(280*w2+220*wl3-340*wl7-160*w26 

-120*w34-60*w221-940))*t4+(36*w2-26*wl3-18*wl7+10*w34-2*w221+6*w442-42)*r2*r3+(74*w2 

+62*wl3+26*wl7-54*w26-2*w34+26*w221+14*w442+286)*r2+(80*w2-90*wl3+10*wl7+40*w26+10*w221 

-130)*r3-470*w2-260*wl3-140*wl7+110*w26-170*w34+20*w221-30*w442-1060; 

=[ cOl, c02, c03, c04, c05, c06, c07, c08, c09, clO, ell, cl2, cl3, cl4, cl5, cl6] ; 
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